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1.  INTRODUCTION 

With  recent  progress  in  integrated  circuit  technology,  MOS  logic 
circuit  has  become  one  of  the  most  important  logic  families  for  digital 
computers.   MOS  logic  circuit  has  a  lot  of  advantages  over  bipolar 
devices  such  as  high  packing  density,  lower  power  consumption  and  simple 
production  process. 

Since  at  least  theoretically  a  MOS  cell  can  realize  an  arbitrary 
negative  function,  several  algorithms  for  designing  logic  circuit  based 
on  MOS  cell's  capability  of  expressing  any  negative  gate  have  been 
developed.   An  algorithm  which  derives  a  two  level  network  of  a  minimum 
number  of  negative  gates  was  first  developed  by  Ibaraki  and  Muroga  [1] 
[2] .   Then  efficient  algorithms  which  can  also  realize  multilevel  network 
with  minimum  number  of  negative  gates  were  developed  by  Liu  [3]  and 
Kasami  et  al  [4] . 

Although  the  algorithms  mentioned  above  guarantee  the  minimality 
of  the  number  of  MOS  cells,  the  network  synthesized  by  these  algorithms 
may  still  have  redundant  connections  and  MOS  FETs.   Recently,  H.  C.  Lai 
[5]  has  modified  those  algorithms  and  introduced  a  new  algorithm  called 
DIMN  (Design  of  Irredundant  MOS  Network)  which  finds  irredundant  MOS 
networks  with  a  minimum  number  of  MOS  cells  for  a  given  set  of  functions. 

In  this  paper,  a  FORTRAN  program  package  DIMN  which  designs  MOS 
logic  networks  based  on  Lai's  algorithm  is  described.   This  program  DIMN 


f 
In  the  case  of  MOS  network,  we  usually  use  the  term  "cell"  in- 
stead of  "gates." 


is  applicable  to  networks  with  multiple  incompletely  specified  output 
functions.   Consequently,  it  covers  almost  all  the  algorithms  introduced 
in  Lai's  paper.   In  other  words,  this  program  has  the  capability  to 
realize  networks  with  a  completely  specified  single  output  function, 
networks  with  completely  specified  multiple  output  functions  and  networks 
with  an  incompletely  specified  single  output  function.   The  only  exception 
is  the  network  with  all  the  output  functions  in  the  output  level.   As 
the  algorithm  for  realizing  this  network  requires  a  slightly  different 
procedure  from  the  procedures  for  realizing  the  networks  mentioned 
above,  it  is  not  implemented  in  the  program  DIMN.   Another  fact  that  has 
to  be  mentioned  here  is  that  the  current  program  DIMN  obtains  only  one 
irredundant  MOS  network  for  any  given  function.   It  is  hoped  that  this 
program  will  be  eventually  modified  to  exhaust  all  the  irredundant  MOS 
networks. 

The  next  chapter,  Chapter  2,  is  allocated  for  the  review  of 
Algorithm  DIMN.   In  Chapter  3,  the  program  DIMN  is  discussed  in  greater 
detail.   Chapter  4  outlines  the  preparation  of  the  input  for  this  program. 
Chapter  5  describes  the  output  of  program  DIMN  and  also  compares  the  net- 
works obtained  by  applying  Algorithm  DIMN  with  the  networks  obtained  by 
Liu's  algorithm.   Finally  the  networks  obtained  by  program  DIMN  for 
several  three  and  four  variable  functions  and  a  complete  listing  of 
FORTRAN  program  DIMN  are  given  in  Appendices  A  and  B,  respectively. 


2.  THEORETICAL  BACKGROUND 

This  section  reviews  Lai's  Algorithm  (Algorithm  DIMN) .   Like  those 
previously  developed  algorithms,  Algorithm  DIMN  repeats  two  phases, 
that  is,  phase  1,  the  derivation  of  a  function  for  each  negative  gate 
(MOS  cell)  and  phase  2,  the  design  of  an  irredundant  MOS  cell  configura- 
tion for  the  function  obtained  in  phase  1.  The  repetition  of  these  two 
phases  leads  to  the  design  of  an  entire  network.   The  difference  between 
Algorithm  DIMN  and  the  previously  developed  algorithms  can  be  seen  in 
how  these  two  steps  are  implemented.   In  the  previously  developed 
algorithms,  the  implementation  is  accomplished  with  a  single  pass  appli- 
cation of  phase  1  followed  by  phase  2.   On  the  contrary,  in  Algorithm 
DIMN,  these  two  phases  are  applied  interactively  to  guarantee  the  irre- 
dundancy  of  the  network. 

The  simplest  case  of  Algorithm  DIMN  i.e.,  the  one  for  obtaining  a 
network  with  a  completely  specified  single  output  function  is  shown  in 
the  following.   This  simplest  version  of  Algorithm  DIMN  can  be  easily 
extended  to  the  cases  for  a  network  with  an  imcompletely  specified 
single  output  function,  for  a  network  with  completely  specified 
multiple  output  functions  and  for  a  network  with  incompletely  specified 
multiple  output  functions.   In  order  to  facilitate  our  discussion, 
some  notations  and  terminologies  are  explained  (see  Lai's  thesis  for 
detail) . 

Our  design  objective  is  to  obtain  a  loopless  network  which  consists 
of  negative  gates  only  for  a  given  function  f.   The  negative  gate  is  a 


gate  which  realizes  a  switching  function  which  can  be  expressed  in  the 

form  of  the  complement  of  a  disjunctive  form  with  non-complemented 

literals  only.  A  generalized  form  of  a  loopless  network  with  R-;  negative 

gates  is  shown  in  Fig.  2.1,  where  x..  , . . .  ,x  denote  N  external  variables 

and  u. , . . . ,u_  denote  the  functions  realized  by  the  negative  gates 
1      Rf 

g-,...,g_   in  the  network. 

NFS  (Rf,f)  =  (u1,...,uR  _v    f): 

In  this  generalized  form  of  a  loopless  network,  the  sequence  of 
functions  u  , . . . ,uR  is  called  a  negative  function  sequence  of  length 

R^  for  a  function  f  and  is  denoted  by  NFS  (R.-,f)  =  (u,  , . . .  ,uD   -  ,  f )  . 
r  r        ±       f~  1 

NFS1  (Rf,f)  =  (ur  . . .  ,u±,  ui+1*,...,uR  _x*,  f): 

This  represents  a  partially  specified  negative  function  sequence 
of  length  R,  and  degree  i  for  a  function  f.  Unlike  the  previously 
mentioned  NFS  (Rf,f),  the  functions  in  this  NFS   (Rf,f)  are  not  com- 
pletely specified.   (In  NFS   (R-,f),  first  i  functions  with  no  *  are 
completely  specified,  but  the  remaining  R-  -  i  -  1  functions  with  * 
are  unspecified.)   A  completion  of  NFS   (R, ,f)  is  a  function  sequence 
obtained  by  completely  specifying  the  unspecified  functions  in 
NFS1  (Rf,f). 

N-cube:   This  is  a  lattice  which  represents  switching  functions 
with  N  external  variables.   In  the  N-cube,  each  vertex  corresponds  to 
an  input  vector  to  the  functions  and  the  vertices  of  the  same  weight 


D — 


Fig.  2.1   Generalized  form  of  a  loopless  network 
consisting  of  Rf  negative  gates. 


(number  of  ones  in  an  input  vector  assigned  to  the  vertex)  are  in  the 
same  level,  placing  vertices  with  more  weight  in  a  higher  level.   Every 
pair  of  vertices  which  corresponds  to  input  vectors  differing  in  only 
one  bit  position  is  connected  by  an  edge.   An  example  is  shown  in 


Fig.  2.2  for  N-3  and  two  functions  (f. 


X..  v  Xq  j    o  X-XaXq  /  • 


Fig.  2.2   3-cube  for  2  functions;  f.  =  x-Vx-, 

t  ry         —   X  -  X^  V  X-XnAn  • 


Algorithm  CMNL:   Algorithm  CMNL  (Conditional  Minimum  Labeling) 
obtains  NFS    (Rt,f)  =  (un , . . . ,u .  ,,  uJt... ,u„    ,  f)  which  is  the 

completion  of  NFS    (R, ,f)  such  that  the  label  assigned  to  each  vertex 
in  the  N-cube  takes  the  minimum  possible  value  of  all  feasible  comple- 
tions of  NFS    (Rf,f).   The  feasible  completion  is  a  completion  such 
that  the  resulting  N-cube  has  no  inverse  edge. 


Algorithm  CMXL:   Algorithm  CMXL  (Conditional  Maximum  Labeling) 
obtains  NFS    (Rf,f)  =  (u, , . . . ,u   , ,  u . , . . . ,uR  _ ,  f)  which  is  the 

completion  of  NFS  (Rf,f)  such  that  the  label  assigned  to  each  vertex 
in  the  N-cube  takes  the  maximum  possible  value  of  all  feasible  comple- 
tions of  NFS1_1(Rf ,f). 

u . :  This  is  the  maximum  permissible  function  which  is  obtained  in 
the  process  which  will  be  described  in  Step  4  of  Algorithm  DIMN. 


Algorithm  DIMN:  Design  of  irredundant  MOS  networks  with  a  minimum 
number  of  MOS  cells  for  a  given  function  f.   Rf  represents  the  minimum 
number  of  MOS  cells. 

Step  1   Let  NFS°(Rf,f)  =  (u*  . . . ,uR   *,  f)  and  set  i=l 

Step  2   Use  algorithm  CMNL  to  obtain  NFS  "  (Rf,f)  =  (u  ,  ...,u.  ,, 

±L"'"\-V   f) 
Step  3   Use  algorithm  CMXL  to  obtain  NFS  "  (Rf,f)  =  (u.,, . . . »u.  ,, 

fii Y1' f) 

Step  4   Obtain  function  u.  by  setting 

u^A)  -  0,  if  u±(A)  =  ui(A)  =  0 

u±(A)  -  1,  if  ^(A)  =  Ui(A)  =  1 

u  (A)  =  *,  if  u^CA)  =  0  and  u±(A)  =  1 

Step  5   Obtain  an  irredundant  MOS  cell  configuration  for  u.  with  re- 
spect to  x. x  ,  u. , . . . ,u .  ,.   Let  u.  denote  the  function 

r       1'     n'   1'    '  i-1       i 

realized  by  this  MOS  cell  (u.  is  now  a  completion  of  u .) . 


Step  6   If  i  =  Rf-1,  design  an  irredundant  MOS  cell  configuration  for 
f  with  respect  to  x..  , . .  .  ,x  ,  u..  , . . .  ,u„  _-  and  terminate  this 

algorithm.   Otherwise  set  i  =  i+1  and  go  to  Step  2. 


It  should  be  noted  that  in  the  above  Algorithm  DIMN,  Steps  2,  3  and  4 
are  deterministic,  in  other  words,  given  a  NFS    (Rf,f),  subsequent 
NFS1"1(Rf  ,f),  NFS1_1(Rf,f)  and  N^Si_:L(Rf  ,f )  are  uniquely  determined. 
On  the  other  hand,  Step  5  in  general  is  non-deterministic,  because 
more  than  one  irredundant  MOS  cell  configuration  may  exist  for  a  given 

V 

Another  fact  which  should  be  noted  is  that  although  Algorithm  DIMN 
requires  at  most  Rf-1  iterations  of  the  loop  which  consists  of 
Step  2,..., Step  6  for  i  =  l,...,Rf-l,  if  NFS    (Rf,f)  becomes  a  com- 
pletely specified  function  with  respect  to  x1,x„,...,x  for  some 
i  <  Rf-1,  the  algorithm  actually  requires  only  i  iterations.   Although, 
even  in  this  case,  Step  5  of  the  algorithm  still  has  to  be  executed 
Rf-i-l  additional  times  in  order  to  obtain  irredundant  MOS  cell  config- 
uration for  u  ,. , . . . ,uD    ,  this  fact  will  help  reducing  the  computation 
1+1       ■p--'- 

time  when  the  algorithm  is  implemented  in  computer  program. 

In  Fig.  2.3,  a  simple  example  is  shown  for  better  understanding  of 
Algorithm  DIMN.  In  this  example,  the  3-cube  with  respect  to  a  function 
f  to  be  realized  is  shown  in  Fig.  2.3(a). 

Step  2  and  Step  3  of  Algorithm  DIMN  obtains  NFS°(3,f)  and  NFS°(3,f) 
as  shown  in  (b)  and  (c)  ,  respectively.   Then  Step  A  compares  u.  and  u 


and  obtains  u  and  NFS°(3,f)  as  shown  in  (d) .   The  NFS1(3,f)1  ■ 
(x  ,  u  *,  f)  in  (e)  is  obtained  by  Step  5.   Step  5  also  obtains  other 
two  irredundant  MOS  cell  configurations  for  u-  which  are  shown  in  (j) 
and  (o) .   After  Step  6,  the  algorithm  returns  to  Step  2  and  obtains 

NFS1(3,f)1,  NFS1(3,f)1,  N^S1(3,f)1  and  NFS2(3,f)1  as  shown  in  (f ) , 

2 
(g),  (h)  and  (i) ,  respectively.   Since  R  =  3,  NFS  (3,f).  is  a  com- 
pletely specified  negative  function  sequence  with  respect  to  x. ,  x„ 
and  x~.   To  finish  the  design,  Step  6  of  the  algorithm  obtains  an 
irredundant  MOS  cell  configuration  for  f .   As  previously  mentioned, 
three  irredundant  MOS  cell  configurations  NFS  (3,f)  ,  NFS  (3,f)2  and 

NFS  (3,f)»  are  obtained  in  (e) ,  (j)   and  (o) ,  respectively  for  given 

a. 

u.  in  (d) .   For  two  other  irredundant  MOS  cell  configurations, 

NFS  (3,f)2,  NFS  (3,f)3,  we  can  continue  applying  Algorithm  DIMN  in  the 
same  way  starting  from  (j)  and  (o) ,  respectively.  The  resulting  irre- 
dundant MOS  networks  are  shown  in  Fig.  2.4  and  Fig.  2.5. 
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(a)   C3(f)  for  f  =  x1x2vx-x 


110 


100 


110 


100 


(b)   NFS°(3,f)  =  (ur  u2,  f) 


110 


100 


(c)   NFSu(3,f)  =  (u1,  u2,  f) 


(d)   NFS°(3,f)  =  (&  ,  u2*,  f) 


110 


100 


110 


100 


(e)   First  NFS1(3,f)]  =  (x^  u2*,  f)    (f)   NFS1(3,f>1  -  (x^  u  2 ,  f) 


Fig.  2.3   Example  for  Algorithm  DIMN. 
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(g)      NFSX(3,f)1  =    (xx,   a2,    f) 


110 


100 


000(110 
(h)      N^S1(3,f)1  =    (5   ,   u2,    f) 


110 


100 


(«    **s*o.th  -  m(3.f)1  -(v2>f)    0>    0-^-™^'Sj  "  <X2,U2*'£) 


110 


(k)      NFSx(3,f)2   =    (x2,   u2,    f) 


(O      NFS    (3,f)2   =    (x2,   u2,    f) 


Fig.    2.3        (Continued) 
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(m)   N^S1(3,f)2  =  (x2,  u2,  f) 


110 


100 


(o)   Third  NFS  (3,f)3=(x3,u2*,f) 
obtained  from  (d) 


110 


100 


110 


100 


(n)   NFS  (3,f)2=NFS(3,f)2=(x2,x1,f) 


110 


100 


(p)   NFS±(3,f)3  =  (x3,  u2,f) 


110 


100 


(q)   NFSi(3,f)3  =  (x3,  u2,  f) 


(r)   N^S1(3,f)3=NFS(3,f)=(x3,f,f) 


Fig.  2.3    (Continued) 
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Fig.  2.4   Irredundant  MOS  network  corresponding  to 
NFS(3,f)'s  in  (i)  and  (n) . 
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Fig.  2.5   Irredundant  MOS  network  corresponding  to 
NFS(3,f)3  in  (r) . 
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3.  PROCEDURE  DIMN 

This  chapter  will  discuss  the  design  procedure  in  FORTRAN  program 
DIMN  for  designing  an  irredundant  MOS  network. 

An  input  to  this  program  is  a  description  of  output  functions  in 
truth  table  form  which  the  resulting  network  has  to  realize.   This 
description  (explained  in  detail  in  Chapter  4)  consists  of  network 
parameters  and  output  function  description.   The  output  of  this  program 
is  the  description  of  realized  irredundant  MOS  networks  (explained  in 
detail  in  Chapter  5) . 

The  entire  DIMN  program  requires  168K  bites  of  core  storage, 
about  62K  being  occupied  by  the  actual  program  instructions  and  106K 
by  the  stored  data  (compiled  by  FORTRAN  G  compiler) . 

For  better  understanding  of  the  procedure  DIMN,  Section  3.1 
describes  internal  data  representation,  that  is,  the  representation  of  a 
labeled  N-cube  (for  the  definition  of  this  terminology  see  Chapter  2) . 
Section  3.2  describes  the  general  organization  of  this  program  and 
Section  3.3  explains  each  subprocedure  in  more  detail.   Although  the 
explanation  of  the  variables  and  arrays  appearing  in  this  program  can 
be  found  in  the  program  listing  in  the  appendix,  Section  3.3  also 
defines  some  of  the  variables  in  order  to  explain  each  flowchart. 

In  the  following  discussion,  N  is  the  number  of  external  vari- 
ables, M  is  the  number  of  output  functions  and  RF  is  the  number  of 
MOS  cells  obtained  by  this  algorithm.   I_I  is  a  pointer  which  indicates 
the  iteration  of  the  program  loop  for  determining  one  MOS  cell 
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configuration.   For  example,  if  TI  is  three,  the  third  MOS  cell  con- 
figuration is  going  to  be  determined. 

3.1  Internal  Data  Representation 

In  order  to  implement  Algorithm  DIMN  in  a  FORTRAN  program,  a 
labeled  N-cube  has  to  be  represented  in  the  computer  memory  effectively, 
This  section  describes  how  the  representation  can  be  accomplished. 

To  implement  the  Algorithm  DIMN  (The  flowchart  of  this  algorithm 
is  shown  in  Fig.  3.1.1.   The  algorithm  shown  in  Fig.  3.1.1  is  slightly 
modified  for  networks  with  incompletely  specified  multiple  output 
functions.)  two  kinds  of  labeled  N-cubes  are  required.   One  is  the 
labeled  N-cube  for  obtaining  MPF  (Maximum  Permissible  Function)  and 
the  other  is  the  labeled  N-cube  for  obtaining  an  irredundant  MOS 
cell  configuration.   The  labeled  N-cube  for  obtaining  MPF  is  accessed 
in  block  2,  block  3  and  block  4  in  Fig.  3.1.1.   In  the  comments  of  the 
program  listing  in  the  appendix,  this  labeled  N-cube  is  simply  referred 
to  as  N-cube  and  each  vertex  is  assigned  a  vertex  number  which  repre- 
sents the  input  vector  in  binary  form. 

As  shown  in  Fig.  3.1.2,  each  vertex  in  a  N-cube  contains  the 

VERTEX  NO. 


LABEL 

DCARE 

MNL 

MXL 

CHAIN 

Fig.  3.1.2   The  vertex  in  the  N-cube. 
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Let  NFS°(Rf,  f 
and  aet  1-1. 


f  *) 


L'se  Algorithm  CMNL  to  obtain 
NFS1"1(R   fj 


fM)-(ul Vl^ %-M'il- 


■..L.). 


Use  Algorithm  CMXL  to  obtain 
NFS1"1(R£.f1 fM)-(Uj 


,f„). 


Obtain  MPF  u   by  6ettlng 


u,(A) 


Obtain  an  lrredundant  MOS  cell  configuration  for 

u"   with  respect  to  x x  ,  u,....,u 

1  1       n    1       1-1 

Let  u   denote  the  function  realized  by  this 


[_ 


Obtain  an  lrredundant  MOS  cell 


Fig.    3.1.1        Algorithm   D1MN, 
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following  five  fields;  LABEL,  DCARE,  MNL,  MXL  and  CHAIN.   LABEL  field, 
together  with  DCARE  field,  stores  the  partially  specified  negative 
function  sequence  (in  Fig.  3.1.1,  this  is  shown  with  the  notation 
NFSi_1(Rf,  f  ,...,fM)).  MNL  and  MXL  fields  store  NFS1~1(Rf,  f  ,...,f  ) 
that  is  the  completion  of  NFS  ~  (Rf,  f  ,...,f  )  by  CMNL,  and 
NFS1"  (R  ,  f  ,...,f  )  that  is  the  completion  of  NFS  "  (Rf,  f  ,...,fM) 
by  CMXL,  respectively.   CHAIN  field  stores  the  link  to  the  next  vertex 
in  the  list  of  the  vertices  of  the  same  weight.   Vertex  0,  that  is, 
the  vertex  with  vertex  number  0  contains  the  first  elements  of  the 
arrays  LABEL,  DCARE,  MNL,  MXL  and  CHAIN;  vertex  1,  that  is,  the  vertex 
with  vertex  number  1  contains  the  second  elements  of  these  arrays  and 
so  on.   In  general,  the  vertex  with  vertex  number  (j-1)  contains  the 
j-th  elements  of  these  arrays  as  shown  in  Fig.  3.1.3. 


VERTEX  (j-1) 


LABEL  (j) 


DCARE  (j) 


MNL  (j) 


MXL  (j) 


CHAIN  (j) 


Fig.  3.1.3   The  vertex  with  vertex  no.  (j-1). 


Fig.  3.1.4  is  an  example  of  N-cube  for  N=3  and  Fig.  3.1.5  is  the 
actual  representation  of  the  3-cube  in  the  computer  memory.   As  can  be 
seen  in  Fig.  3.1.4,  the  first  and  the  last  vertices  in  the  lists  of  the 
vertices  with  the  same  weight  are  pointed  by  pointers  STARTL  and  ENDL, 
respectively.   For  example,  vertex  3,  the  first  vertex  in  the  list  of 
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LABEL 


DCARE 


MNL 


MXL 


CHAIN 


(1) 
(2) 
(3) 
(4) 
(5) 
(6) 
(7) 
(8) 


(1) 
(2) 
(3) 
(4) 


STARTL 


ENDL 


Fig.  3.1.5   Internal  representation  of  the  3-cube  in 
Fig.  3.1.4. 


the  vertices  with  weight  2  is  pointed  by  pointer  STARTL(3)  and 
vertex  6,  the  last  vertex  in  the  same  list  is  pointed  by  pointer 
ENDL(3) .   In  general,  the  first  and  the  last  vertices  in  the  list  of 
the  vertices  with  weight  K  is  pointed  by  pointers  STARTL  (k+1)  and 
ENDL  (k+1),  respectively. 

The  labeled  N-cube  for  obtaining  an  irredundant  MOS  cell  configura- 
tion is  accessed  in  block  5  in  Fig.  3.1.1.   In  the  comments  of  the 
program  listing  in  Appendix  B,  this  labeled  N-cube  is  referred  to  as 
the  Large-cube  because  the  dimension  N  of  this  N-cube,  starting  from 
the  initial  value  N,  is  increased  by  one  every  time  the  program  loop 
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in  Fig.  3.1.1  is  executed.   The  above  discussion  implies  that  this 
Large-cube  has  to  be  newly  constructed  every  time  block  5  in 
Fig.  3.1.1  is  executed.   On  the  other  hand,  the  previously  mentioned 
N-cube,  once  being  constructed  at  the  beginning  of  the  program,  will 
never  be  executed  (the  contents  of  the  N-cube  is  changed  every  time 
the  program  loop  in  Fig.  3.1.1  is  executed). 

In  a  similar  way  to  the  N-cube,  each  vertex  in  the  Large-cube 
which  consists  of  FUNC  and  LINK  fields  is  assigned  a  vertex  number 
which  is  the  decimal  representation  of  the  corresponding  input  vector. 
The  LINK  field  corresponds  to  the  CHAIN  field  of  the  N-cube  and  points 
to  the  next  vertex  with  the  same  weight.  MPF  (Maximum  Permissible 
Function)  is  stored  in  the  FUNC  field  after  MPF  is  obtained  at  the 
end  of  block  4  in  Fig.  3.1.1. 

At  the  i-th  iteration  of  the  program  loop  in  Fig.  3.1.1,  the 
dimension  of  this  Large-cube  will  be  N+i-1  and  every  time  the  program 
loop  is  executed,  the  dimension  will  be  increased  by  one.   Therefore, 
at  the  last  iteration  of  this  program  loop  (when  i  ■  Rf ) ,  the  dimension 

of  this  Large-cube  will  be  N  +  Rf  -  1.   This  means  that  when  N=10  and 

14 
Rf  =  5,  the  memory  space  for  storing  2   vertices  is  required  in  order 

to  store  this  Large-cube.   (The  restriction  on  problem  size  is  due  to 

this  memory  size.   This  will  be  discussed  in  Section  4.2.) 

Fig.  3.1.6  shows  an  example  for  3-dimensional  Large-cube  (this 

means  N  +  i  -  1  =  3)  and  Fig.  3.1.7  is  the  internal  representation  of 

the  3-dimensional  Large-cube  shown  in  Fig.  3.1.6. 


STARTF(4)    111 


ENDF(A) 
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ENDF 


Fig.  3.1.7   Internal  representation  of  3 -dimensional 
Large-cube  in  Fig.  3.1.6. 
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3.2  General  Organization  of  Program  DIMN 

This  section  shows  the  general  organization  of  program  DIMN  and 
outlines  each  subprogram.   The  function  of  each  subprogram  is  discussed 
in  detail  in  the  following  section. 

The  general  organization  of  program  DIMN  is  shown  in  Fig.  3.2.1. 
As  can  be  seen  in  this  diagram,  this  program  consists  of  subroutine 
MAIN  and  the  following  subroutines:  INPUT,  NCUBE,  CMNL,  CMXL,  MPF, 
IMC,  ASFUN1,  ASFUN2,  INCMNT,  DSCAN,  ASIGN1,  DCRMNT,  CMPR  and  PRINT. 
In  Fig.  3.2.1,  an  arrow  from  block  i  to  block  j  represents  the  fact 
that  the  subroutine  represented  by  block  i  calls  the  subroutine 
represented  by  block  j . 

Subroutine  INPUT:   This  subroutine  reads  into  an  N-cube  input 
data,  that  are  problem  parameters  and  given  output  functions.   Input 
data  setup  is  described  in  detail  in  Chapter  A. 

Subroutine  NCUBE:   This  subroutine  constructs  an  N-cube  according 

to  the  number  of  external  variables.   This  means  if  the  number  of 

N 
external  variables  is  N,  this  N-cube  will  contain  2  vertices. 

Subroutine  NCUBE  will  examine  the  input  vector  in  binary  form  which 
is  assigned  to  each  vertex  in  the  N-cube  one  by  one  and  link  the 
vertices  with  the  same  weight  together. 

Subroutine  CMNL:   This  subroutine  implements  Conditional  Minimum 
Labeling  in  Algorithm  DIMN  based  on  the  output  function  values  read 
into  the  N-cube.   As  shown  in  Fig.  3.2.1,  this  subroutine  calls  sub- 
routines INCMNT  and  DSCAN.   RF,  the  minimum  number  of  MOS  cells  which 


23 


CM 

tn 

CO 

<: 

55 

i 

E 

CO 

<3 

H 

PC! 
Ph 

P4 

" 

Pn 

2 

g 

U 

53 

O 

M 

i-H 

S3 

H 

O 

-..g 

M 
CO 

<d 

l-J 

»! 

o 

CO 

u 

o 

w 

H 

W 

^U 

53 

55 
H 

H 

P 

^ 

M 

CM 

cn 

•H 


24 


is  required  to  realize  the  given  function(s)  is  determined  at  the  first 
execution  of  this  subroutine  CMNL. 

Subroutine  CMXL;   This  subroutine  implements  Conditional  Maximum 
Labeling  in  Algorithm  DIMN  based  on  the  output  function  values  read 
into  the  N-cube.  This  subroutine  calls  subroutine  DSCAN,  ASIGN1  and 
DCRMNT. 

Subroutine  MPF:   This  subroutine  obtains  the  maximum  permissible 
function  from  the  results  obtained  in  MNL  and  MXL  fields  in  the  N-cube 
by  subroutines  CMNL  and  CMXL,  respectively.   Then  this  subroutine 
stores  the  resulting  MPF  into  the  FUNC  field  of  the  Large-cube. 

Subroutine  IMC:   This  subroutine  obtains  an  irredundant  MOS  cell 
configuration  from  the  maximum  permissible  function  in  the  Large-cube 
obtained  by  subroutine  MPF.   This  subroutine  is  constructed  with  the 
following  six  steps.   This  subroutine  calls  subroutines  COMPR  and  PRINT. 

Step  1:   In  this  step,  the  Large-cube  is  constructed  in  the  same  way 

as  the  N-cube  is  constructed  in  subroutine  NCUBE.   The  dimen- 
sion of  this  Large-cube  is  determined  by  the  number  N  +  RF  -  1. 

Step  2:   This  step  assigns  "0"  or  "1"  to  each  vertex  X  which  has  don't 
care  in  the  Large-cube  such  that  there  exists  no  vertex  Y 
satisfying  the  following  condition: 

Input  vector  Y  >  Input  vector  X  and  the  function  value 
assigned  to  the  vertex  Y  is  "1." 

Step  3:   This  step  obtains  the  set  of  minimum  vectors. 

Step  4:   This  step  obtains  a  subset  of  the  set  of  minimum  vectors  which 
covers  all  the  vertices  which  originally  have  the  value  "0." 
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Step  5:  This  step  obtains  an  irredundant  subset  from  the  subset  obtained 
in  Step  4.  The  irredundant  subset  obtained  in  this  step  will 
be  used  for  the  realization  of  an  irredundant  MOS  cell  con- 
figuration.  Since  the  finding  of  all  irredundant  covers 
(subsets)  is  essentially  the  covering  problem,  it  will  be 
very  time-consuming  as  the  problem  size  becomes  larger.   So 
in  this  subroutine,  only  one  irredundant  subset  is  obtained. 

Step  6:   This  step  stores  values  of  the  function  realized  by  the  MOS 

cell  obtained  in  Step  5  in  the  LABEL  field  in  the  N-cube.   We 
have  to  be  careful  that  the  function  value  realized  at  each 
vertex  with  don't  care  by  the  irredundant  MOS  cell  may  be 
different  from  the  function  value  assigned  to  the  corresponding 
vertex  in  the  Large-cube  in  Step  2. 

This  step  also  contains  an  error  checking  routine.   If  the 
function  values  realized  by  the  irredundant  MOS  cell  obtained 
in  Step  5  are  different  from  those  already  stored  in  the  LABEL 
field  of  the  N-cube,  an  error  message  and  contents  of  the 
array  LABEL  is  printed. 

Subroutine  ASFUN1:   Whenever  certain  conditions  which  will  be  dis- 
cussed in  detail  in  Section  3.3,  are  met,  Algorithm  DIMN  does  not  need 
to  obtain  the  maximum  permissible  function.   This  means  that  the  program 
loop  which  consists  of  CMNL,  CMXL  and  MPF  need  not  be  implemented  and  by 
detecting  such  conditions,  we  can  save  computation  time.   This  sub- 
routine is  called  under  such  circumstances.   Then  the  function  values 
stored  in  the  LABEL  field  in  the  N-cube  are  directly  stored  by  this 
subroutine  in  the  FUNC  field  in  the  Large-cube. 
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Subroutine  ASFUN2:   This  subroutine  is  called  under  a  similar  con- 
dition to  that  subroutine  ASFUN1  is  called.   The  values  in  the  MNL 
field  in  the  N-cube  are  stored  by  this  subroutine  in  FUNC  field  in 
the  Large-cube. 

Subroutine  MAIN:   This  subroutine  calls  the  subroutines  described 
above,  whenever  necessary,  and  implements  Algorithm  DIMN. 

3.3  Description  About  Subroutines 

(1)   Subroutine  MAIN  (Fig.  3.3.1) 

In  block  1,  parameters  N  and  M  are  read  in  and  if  these  parameters 
do  not  satisfy  the  restriction  on  problem  size  (see  Section  A. 2  for 
detailed  discussion) ,  an  error  message  is  printed  out  and  the  program 
execution  is  terminated.   This  block  also  tests  EOF  (end  of  file)  con- 
dition and  if  all  data  have  been  exhausted,  the  program  execution  is 
terminated . 

In  block  2,  the  LABEL  and  DCARE  fields  in  the  N-cube  are  initial- 
ized to  zero. 

In  block  3,  the  subroutine  INPUT  is  called  and  values  on  output 
function  cards  are  read  into  the  LABEL  field  in  the  N-cube. 

In  block  4,  the  N-cube  is  constructed  based  on  the  parameter 
value  N. 

In  block  5,  II ,  the  pointer  which  indicates  the  number  of  iterations 
of  the  program  loop  is  initialized  to  1.   One  MOS  cell  configuration  is 
determined  every  time  the  loop  is  executed. 
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CALL  ASFIN1 

»                 , 
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CALL   IMC 

«                , 
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II 

CALL    IMC 

Fig.  3.3.1   Flowchart  of  subroutine  MAIN. 
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Block  6  implements  conditional  minimum  labeling  (CMNL)  and  at  the 
first  execution  of  this  block  the  value  of  RF  (the  minimum  number  of 
MOS  cells)  is  determined. 

In  block  7 ,  the  obtained  RF  value  is  examined  and  if  it  is  1,  the 
maximum  permissible  function  (MPF)  need  not  be  obtained;  the  MOS  cell 
configuration  can  be  obtained  immediately  by  calling  ASFUN1  and  IMC. 
If  R  ^  1,  then  subroutine  CMXL  is  called  to  find  the  conditional 
maximum  labeling. 

In  block  9,  after  executing  subroutine  CMXL,  the  labels  assigned 
to  the  N-cube  by  CMNL  and  CMXL  are  compared  and  if  the  labels  assigned 
by  the  two  different  ways  take  the  same  value  at  every  vertex  in  N- 
cube,  there  exists  unique  minimum  negative  gate  network  and  block  9 
is  immediately  followed  by  the  loop  for  obtaining  an  irredundant  MOS 
cell  configuration.   Under  this  circumstance,  the  sequence  for  obtain- 
ing maximum  permissible  function,  that  is,  the  sequence  of  subroutines 
CMNL,  CMXL  and  MPF  is  skipped. 

In  block  9,  if  the  labels  assigned  to  the  N-cube  take  different 
values  at  some  vertices,  this  block  is  followed  by  block  10  (subroutine 
MPF)  and  block  11  (subroutine  IMC) .   After  the  implementation  of  block 
11,  the  configuration  of  MOS  cell  indicated  by  pointer  LI  is  deter- 
mined (the  II-th  MOS  cell  is  obtained) . 

In  block  12,  pointer  LI  is  increased  by  one  and  in  block  14,  if 
II  is  not  equal  to  RF,  in  other  words,  if  we  are  not  going  to  deter- 
mine the  last  MOS  cell  configuration,  the  program  control  is  returned 
to  block  6.   On  the  contrary,  if  II  is  equal  to  RF,  that  is,  if  the 
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last  MOS  cell  configuration  is  going  to  be  determined,  the  procedure 
for  determining  MOS  cell  configuration  is  immediately  followed  (the 
procedure  for  obtaining  MPF  is  skipped) .  After  the  implementation  of 
block  16,  an  irredundant  MOS  network  which  realizes  the  given  output 
functions  is  obtained. 

In  block  17 ,  the  statistics  about  the  obtained  network  are  printed 
out  and  the  program  control  is  transferred  to  block  1.  The  contents 
of  the  statistics  are  described  in  detail  in  Chapter  5. 

In  block  1,  if  problem  cards  are  not  exhausted,  a  new  problem  is 
read  in  and  the  entire  process  is  repeated.   If  problem  cards  have 
already  been  exhausted,  the  program  is  terminated. 

In  block  13,  the  logical  variable  CS  is  examined  (this  variable 
is  set  in  subroutine  INPUT  if  all  the  output  functions  are  completely 
specified)  and  if  CS  is  true,  this  block  is  followed  by  block  22.   In 
the  case  of  multiple  output  functions,  the  maximum  permissible  function 
need  not  be  obtained  for  the  output  function  which  is  completely 
specified.  Therefore,  if  one  of  the  MOS  cells  which  realize  the  com- 
pletely specified  given  output  function  is  being  determined,  block  22 
is  immediately  followed  by  the  loop  for  determining  MOS  cell  configura- 
tion and  the  process  for  obtaining  the  maximum  permissible  function  is 
skipped. 

(2)   Subroutine  INPUT 

At  the  beginning  of  this  section,  the  way  how  the  output  functions 
are  stored  in  arrays  LABEL  and  DCARE  will  be  discussed  and  then  flow- 
chart will  be  explained. 
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As  described  in  Section  3.1,  a  partially  specified  negative  func- 
tion sequence  is  stored  in  the  LABEL  and  DCARE  fields  of  the  N-cube. 
The  internal  representation  of  the  LABEL  and  DCARE  fields  in  the  N-cube 
which  has  N  external  variables  and  M  output  functions  is  shown  in 
Fig.  3.3.2.   As  shown  in  this  figure,  the  labels  assigned  to  vertices 
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u.  -  represents  the  function  realized  at  the  output  of  the  i-th 
MOS  cell. 

f   -  represents  output  function  f  . 


Fig.  3.3.2   The  internal  representation  of  the  LABEL  and 
DCARE  fields  in  the  N-cube. 


in  the  N-cube  are  stored  in  bit  pattern.   The  column  labeled  with  u. 
stores  the  function  realized  at  the  output  of  the  i-th  MOS  cell  and 
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the  column  labeled  with  f .  stores  the  given  output  function  f . .  As  will 
be  discussed  in  Chapter  4,  since  the  output  functions  are  supplied  in 
truth  table  form,  output  function  1  is  stored  in  column  f  ,  output 
function  2  is  stored  in  column  f?  and  so  forth.   In  general,  output 
function  i  is  stored  in  column  f . . 

The  function  value  of  the  i-th  function  is  stored  in  the  LABEL  and 
DCARE  fields  of  vertex  (j-1)  in  the  following  way.   (As  described  in 
Section  3.1,  the  LABEL  and  DCARE  fields  of  vertex  (j-1)  are  stored  in 
LABEL  (j)  and  DCARE  (j),  respectively.) 

If  a  function  value  is  zero,  zero  is  stored  by  assigning  zero  to 
column  f  of  both  LABEL  (j)  and  DCARE  (j).   (Since  arrays  LABEL  and 
DCARE  have  already  been  initialized  to  zero  in  block  2  of  subroutine 
MAIN,  no  action  need  be  taken.)   If  a  function  value  is  one,  one  is 
stored  by  assigning  one  to  column  f..  of  LABEL  (j)  and  zero  to  column 
f.  of  DCARE  (j).   This  is  implemented  by  adding  2M_i  to  LABEL  (j). 
Finally,  if  a  function  value  to  be  stored  is  don't  care,  this  is  done 
by  assigning  zero  to  column  f..  of  LABEL  (j)  and  one  to  column  f.  of 
DCARE  (j).   This  is  implemented  by  adding  2M_i  to  DCARE  (j).   The 
example  in  Fig.  3.3.3  shows  how  the  output  functions  in  the  right 
hand  side  are  stored  in  arrays  LABEL  and  DCARE  for  N=3  and  M=2.   In 
the  truth  table  of  this  figure,  *  represents  don't  care  condition. 

The  following  is  the  definition  of  each  variable  which  appears 
in  subroutine  INPUT. 

NCARD:   This  stores  the  number  of  input  data  cards  required  to 
describe  one  output  function. 
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Fig.  3.3.3   Example:   how  output  functions  are  stored 
in  arrays  LABEL  and  DCARE. 


I_:   This  is  an  indicator  of  output  function  number.   If  problem 
requires  M  output  functions,  I  changes  from  one  to  M. 

J^:   This  is  an  indicator  of  card  number  within  one  output  function, 
J  changes  from  one  to  NCARD. 

K:   This  is  a  column  indicator  in  one  card.   Therefore,  K  changes 
from  one  to  80. 

VTEX:   This  is  a  pointer  to  a  vertex  in  the  N-cube  in  which  an 
output  function  value  is  to  be  stored.   When  VTEX  =  i,  this  points 
to  vertex  (i-1) . 

CHAR  (80) :   This  array  is  a  character  buffer  for  input  data. 
This  implies  that  one  output  function  card  is  read  in  at  a  time. 

CS :   This  is  a  logical  variable  which  is  set  when  every  output 
function  is  completely  specified. 
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Fig.  3.3.4  shows  the  flowchart  of  subroutine  INPUT.   In  block  1, 
CS  is  set.   If  don't  care  is  found  in  output  functions,  it  will  be  reset 
in  block  9.  This  means  that  after  executing  subroutine  INPUT,  CS  is  set 

if  and  only  if  every  output  function  is  completely  specified. 

N 
In  block  2,  NCARD  is  obtained  in  the  following  way.   If  2  is 

devisable  by  80,  then  NCARD  =  2N/80,  otherwise  NCARD  =  2N/80+l,  where 

/  represents  integer  division  and  the  result  of  2  /80  is  equal  to 

L-80  -I  * 

In  block  7,  a  function  value  which  is  stored  in  CHAR(K)  is  examined 
and  the  function  value  is  stored  in  the  N-cube  in  the  way  described  in 
the  beginning  of  this  section.   In  the  input  format,  a  blank  means  the 
termination  of  the  description  of  one  output  function.  Therefore,  if  a 
character  in  CHAR(K)  is  a  blank,  this  is  interpreted  as  the  termination 
of  the  description  of  one  output  function  and,  in  block  15,  an  error  is 
checked.   In  a  normal  termination  of  the  description  of  one  output 
function,  at  the  instant  when  a  blank  character  has  been  received,  J 
and  VTEX  have  to  take  the  following  value; 

J  =  NCARD,   VTEX  =  2N  +  1 

If  an  error  occurs,  an  error  message  is  printed  and  the  control 
returns  to  block  1  in  subroutine  MAIN.  This  means  that  this  problem 
has  been  skipped.  In  block  7,  receiving  characters  other  than  zero, 
one,  *  and  blank  is  also  interpreted  as  an  error. 
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3.3.4   Flowchart  of  subroutine  INPUT. 
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(3)  Subroutine  NCUBE 

Flowchart  is  shown  in  Fig.  3.3.5.  The  following  is  the  definitions 
of  variables  which  appear  in  this  subroutine. 

I_:  This  stores  an  input  vector  in  binary  form. 

X:  This  is  a  variable  in  which  each  input  vector  is  examined  bit 
by  bit  and  the  weight  of  the  input  vector  is  determined. 

W:  This  is  a  variable  in  which  the  weight  of  the  input  vector 
is  stored. 

STARTL,  ENDL  and  other  arrays  are  described  in  Section  3.1. 

N 
When  the  number  of  external  variables  is  N,  there  exists  2  input 

N 
vectors  (ranging  from  0  to  2  -1  in  binary  form).  Therefore,  the  pro- 
cedure which  determines  the  weight  of  an  input  vector  and  links  the 

N 
input  vectors  with  the  same  weight  is  repeated  for  1=0  to  2  -1.   In 

the  loop  which  consists  of  blocks  4,  5,  6  and  7,  each  input  vector  is 
examined  bit  by  bit  and  in  blocks  8,  9  and  10,  the  input  vectors  with 
the  same  weight  are  linked  together.   For  this  chain  constructing  pur- 
pose, STARTL  has  to  be  initialized  to  zero  in  block  1.  An  example  of 
the  N-cube  constructed  in  this  procedure  is  shown  in  Fig.  3.1.5  for 
N=3. 

(4)  Subroutine  CMNL 

Flowchart  is  shown  in  Fig.  3.3.6.  The  following  is  the  definition 
of  variables  which  appear  in  this  subroutine. 

USPFY:   As  previously  mentioned,  array  LABEL  stores  a  partially 
specified  negative  function  sequence.  USPFY  stores  the  number  of  un- 
specified bits  in  array  LABEL.   The  relation  between  USPFY,  RF,  II  and 
M  is  shown  in  Fig.  3.3.7. 
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INITIALIZE    START1. 


I    •    0 


X   -    I 
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is. 


CHAIN    (ENDL(W+1)>    *■  I  +  1 

END(W  +    1)    •    1+1 


Fig-  3.3.5   Flowchart  of  subroutine  NCUBE. 
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Fig.  3.3.6   Flowchart  of  subroutine  CMNL. 
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Fig.  3.3.7   An  element  of  array  LABEL, 


W:   This  stores  an  actual  weight  plus  one. 

PTRB:   This  is  a  pointer  to  vertex  B  to  which  the  minimum  possible 
value  is  to  be  assigned. 

PTRA:   This  is  a  pointer  to  vertex  A  which  is  connected  to  vertex 
B  by  the  edge  from  vertex  A  to  vertex  B.   The  vertex  A  has  a  greater 
weight  than  vertex  B  by  one. 

LBX:  This  is  a  variable  in  which  LABEL  pointed  by  pointer  B  is 
increased. 

LAX:   This  is  a  variable  which  stores  the  minimum  label  which  has 
already  been  assigned  to  vertex  A. 

LBY:  This  is  a  variable  in  which  DCARE  pointed  by  pointer  B  is 
examined  in  subroutine  DSCAN.  The  number  of  don't  care  bits  and . the 
weight  assigned  to  each  don't  bit  are  determined. 

NDCARE :   This  is  a  variable  which  stores  the  number  of  don't  care 
bits. 

BWEIT:   This  is  an  array  which  stores  the  weight  assigned  to  each 
don't  care  bit. 

PTRBX:   This  is  a  variable  in  which  the  input  vector  in  binary 
form  assigned  to  vertex  B  is  examined  bit  by  bit.   This  variable, 
together  with  variable  SHTFT,  determine  PTRA  by  calculation. 


39 


INC:  This  is  a  counter  for  the  increment  of  LBX.  This  is 
described  in  detail  in  subroutine  INCMNT. 

In  block  1  (Fig.  3.3.6),  the  minimum  possible  label  is  assigned 
to  a  vertex  with  weight  N.   Since  LABEL  is  initialized  to  zero  in  block  2 
in  subroutine  MAIN  (Fig.  3.3.1),  the  unspecified  bits  in  LABEL  (STARTL 
(N+l))  have  already  been  assigned  zero,  the  minimum  possible  value.  As 
shown  in  Fig.  3.3.7,  the  number  of  unspecified  bits  is  obtained  by 
calculating  RF  -  II  +  1.  As  RF  is  obtained  after  the  first  implementa- 
tion of  subroutine  CMNL,  the  value  of  RF  has  to  be  specified  in  subroutine 
MAIN  before  CMNL  is  called  for  the  first  time. 

After  assigning  a  label  to  the  vertex  with  weight  N,  a  minimum 
possible  label  is  assigned  to  every  other  vertex  in  the  following  way. 
In  blocks  4,  5,  . . . ,  16,  each  vertex  pointed  by  PTRB  (vertex  B)  is 
assigned  a  minimum  possible  label.   In  general,  as  shown  in  Fig.  3.3.8, 

VERTEX  (PTRA-1) 
(      1 (       J [         I  PTRA 


(    J      PTRB 
VERTEX (PTRB-1) 

Fig.  3.3.8   Relation  between  vertices  A  and  B  in  CMNL. 

there  exist  several  vertices  which  are  directly  connected  to  vertex  B 
by  the  edges  from  these  vertices  to  vertex  B.  Minimum  labels  which 
have  already  been  assigned  to  these  vertices  are  compared  with  LBX 
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one  by  one  (block  11,  LABEL  (PTRB)  is  assigned  to  LBX  in  block  A)  and 
LBX  is  increased  by  one  (block  13)  until  the  condition  that  no  one  of 
the  labels  assigned  to  these  vertices  is  larger  than  the  label  assigned 
to  vertex  B  is  satisfied.   In  block  9,  PTRA  is  calculated  by  PTRA  = 
PTRB  +  SHIFT.   As  vertices  with  the  same  weight,  say  (w-1) ,  are  linked 
by  CHAIN  field,  starting  from  the  vertex  pointed  by  STARTL(W),  all  the 
vertices  within  the  linked  list  can  be  exhausted  by  utilizing  this 
scheme  (blocks  3,  17  and  18) . 

Starting  from  the  vertices  with  weight  (N-l) ,  the  above  procedure 
is  applied  to  every  group  of  vertices  with  the  same  weight  until  a 
label  is  assigned  to  the  vertex  with  weight  zero  (blocks  2,  19,  20) . 
At  the  first  implementation  of  this  subroutine,  RF  value  is  obtained 
(blocks  21,  22). 

(5)   Subroutine  DSCAN 

Flowchart  is  shown  in  Fig.  3.3.9.   This  subroutine  has  four 
parameters;  two  of  them  are  input  parameters  (LBY,  USPFY)  and  the 
others  are  output  parameters  (NDCARE,  BWEIT) .   It  is  called  by  CMNL 
and  CMXL. 

It  determines  the  number  of  don't  care  bits  in  the  function  part 
of  the  label  assigned  to  vertex  B  and  the  weight  assigned  to  these 
don't  care  bits.   This  is  accomplished  by  examining  LBY  (don't  care 
field  of  vertex  B)  bit  by  bit.   If  "1"  is  encountered,  NDCARE  is  in- 
creased by  one  and  the  corresponding  bit  weight  is  stored  in  BWEIT 
(NDCARE).   This  bit  weight  is  represented  in  variable  SHIFT  (this 
stores  the  value  2   ) .   The  above  procedure  Is  continued  until  all  the 
unspecified  bits  are  examined. 
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KDCARI  •  0 
SHIFT  «■  1 


I  -  1 


NDCARE  -  NDCARE  +  1 
BWEIT(!TOCARE) -SHIFT 


I  -  I  +  1 

SHIFT  -  2  *  SHIFT 


YES 


LBY  *  LBY/ 2 


YES 


END 


Fig.  3.3.9   Flowchart  of  subroutine  DSCAN, 
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An  example  is  shown  in  Fig.  3.3.10  for  a  label  assigned  to  vertex 

B.   The  label  in  this  example  contains  2  don't  care  bits  and  the  weight 

0      2 
of  these  bits  are  2  and  2  . 


(6)   Subroutine  INCMNT 

This  subroutine  has  five  parameters;  three  of  them  (NDCARE,  BWEIT, 
M)  are  input  parameters  and  the  others  (INC,  LBX)  are  input-output 
parameters.   The  flowchart  is  shown  in  Fig.  3.3.11. 

Although  INC  is  a  counter  for  incrementing  LBX  (the  label  field 

of  vertex  B) ,  each  bit  in  counter  INC  has  a  different  weight  from  an 

ordinary  binary  counter.   Bit  weights  of  the  least  significant  half 

which  correspond  to  the  bit  weights  of  don't  care  bits  in  LBX  can  be 

obtained  in  array  BWEIT  by  subroutine  DSCAN  (the  number  of  these  bits 

is  stored  in  NDCARE) .   Weights  of  other  unspecified  bits  can  be  obtained 

by  the  following  formula;  the  bit  weight  of  the  J-th  least  significant 

bit  is 

2M  +  J  -  NDCARE  -  1 

where  J  is  the  index  which  indicates  the  number  of  program  loop  itera- 
tions in  Fig.  3.3.11.   An  example  of  the  bit  weight  assignment  to 
counter  INC  is  shown  in  Fig.  3.3.10  for  a  label  of  vertex  B. 

In  block  1  (Fig.  3.3.11),  counter  INC  which  has  already  been 
initialized  to  zero  in  block  4  in  Fig.  3.3.6  (CMNL)  is  increased  by 
one. 

In  block  3,  variable  SHIFT  and  index  J  are  initialized  so  that 
SHIFT  stores  the  bit  weight  of  the  J-th  least  significant  bit.   In 
the  loop  which  consists  of  blocks  4, . . . ,10,  the  contents  of  counter 
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.  (vertex  B) 
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,M+2 
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2M    22    20 

The  weights  assigned  to  bits 


Fig.  3.3.10   Example:   how  don't  care  bits  in  the  LABEL 
field  is  treated. 


START 
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i      V 

INC  -  INC  +  1 

2         V 

INCX  *  INC 

3         V 

SHIFT  +   2M-NDCARE 

J  '-   1 

•fc 

-  0 

INCX 

LBX  «.  LBX  +  BWEIT(J) 


RETURN 


LBX  «■  LBX  +  SHIFT 


J  «-  J  +  1 
SHIFT  *■  2  *  SHIFT 


10 


INCX  «-  INCX/2 


Fig.  3.3.11   Flowchart  of  subroutine  INCMNT. 
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INC  are  examined  bit  by  bit  and  if  "1"  is  encountered,  the  weight 
assigned  to  this  bit  of  1  is  added  to  LBX. 

In  block  6,  the  value  of  index  J  is  examined  and  if  J  is  not 
greater  than  NDCARE,  the  contents  of  BWEIT(J)  is  added  to  LBX,  and 
otherwise  the  contents  of  the  variable  SHIFT  is  added. 

After  executing  this  subroutine,  LBX  is  increased  to  the  next 
possible  state  and  then  compared  with  LAX  in  block  11  of  subroutine 
CMNL  (Fig.  3.3.6).   It  is  to  be  noted  that  LBX  is  not  always  increased 
by  one  because  it  may  contain  already  specified  bits.  As  shown  in 
block  12  (Fig.  3.3.6),  prior  to  calling  this  subroutine,  LBX  has  to 
be  initialized  to  the  value  stored  in  LABEL  (PTRB) . 

(7)   Subroutine  CMXL 

This  subroutine  is  executed  in  almost  the  same  way  as  the  sub- 
routine CMNL  is.  The  only  difference  is  that  after  l's  are  assigned 
to  the  unspecified  bits  of  label  field  of  vertex  B  by  calling  sub- 
routine ASIGN1,  the  label  is  decreased  instead  of  being  increased  in 
subroutine  CMNL. 

The  flowchart  is  shown  in  Fig.  3.3.12  and  the  following  is  the 
definitions  of  variables  used  in  this  subroutine.  Variables  with  the 
same  definition  as  those  in  subroutine  CMNL  are  omitted. 

PTRB:   This  is  a  pointer  to  vertex  B  to  which  the  maximum  possible 
value  is  to  be  assigned. 

PTRA:   This  is  a  pointer  to  vertex  A  which  is  connected  to  vertex 
B  by  the  edge  from  vertex  B  to  vertex  A.   The  vertex  A  has  a  smaller 
weight  than  vertex  B  by  one.   The  relation  between  vertex  A  and  vertex 
B  is  shown  in  Fig.  3.3.13. 


LBX  -  LABEL(STARTLU)) 
LBY  »  DCARE(STARTLd)) 
U3PFY   -    RP   -    II   +   1 


CALL   ASIGN1 


HXL(STARTLU))   *  LBX 


PTRB   *   STARTL(W) 


LBY   *  DCARE(PTRB) 


LUX   *  LABEL(PTRB) 
LBY   .   DCARE(PTRB) 


CALL   ASIGN1 


t~ 


LBXX   -    LBX 

INC   -   0 

i  *  l,  shift  •  i 


PTRB   -    CHAIN(PTRB)  ^ 


rrux  •  (PTU-i) 
/  SHirr 


PTRA  *  PTRB   -   SHIFT 


LAX  -  MXL    (PTRA) 


1 5^-^''^~"~~-».^ 

LBXX  <  LAX 

16                 y 

LBB  -  LBX 

17                    \ 

' 

-LS_ 


1*1  +  1 

SHIFT   •  2   •   SHIFT 


MXL(PTRB)    -    LBXX 
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Fig.  3.3.12    Flowchart  of  subroutine  CMXL. 
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VERTEX  (PTRB-1) 
PTRB 


PTRA 
VERTEX  (PTRA-1) 

Fig.  3.3.13   Relation  between  vertices  A  and  B  in  CMXL, 


LBX:  This  is  a  variable  where  LABEL  (PTRB)  is  stored  after  l's 
are  assigned  to  the  unspecified  bits  in  LABEL  (PTRB) . 

LBXX:   This  is  a  variable  in  which  contents  of  LBX  is  determined. 

LAX:   This  is  a  variable  which  stores  the  maximum  label  which  has 
already  been  assigned  to  vertex  A. 

INC:  This  is  a  counter  for  decreasing  LBXX.  The  contents  of  INC 
is  subtracted  from  LBXX  in  subroutine  DCRMNT. 

In  blocks  1,  2  and  3  (Fig.  3.3.12),  the  maximum  possible  label  is 
assigned  to  the  vertex  with  weight  zero.   Starting  from  the  vertices 
with  weight  one  (block  2) ,  a  maximum  possible  label  is  assigned  to 
every  vertex  in  the  N-cube  in  a  similar  manner  as  in  subroutine  CMNL. 
In  block  13,  PTRA  is  obtained  by  subtracting  SHIFT  from  PTRB. 

(8)   Subroutine  ASIGN1 

This  subroutine  has  four  parameters;  three  of  them  (USPFY,  M  and 
LBY)  are  input  parameters  and  the  remaining  one  (LBX)  is  an  input- 
output  parameter.   One  is  assigned  to  every  unspecified  bit  of  LBX 
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including  unspecified  don't  care  bits  in  the  values  of  output  functions 
assigned  to  vertex  B.   This  is  accomplished  by  examining  input  parameter 
LBY. 

In  Fig.  3.3.14,  I   is  a  variable  which  points  to  a  bit  in  the  label 
field  of  vertex  B.   Starting  from  the  least  significant  bit  of  the  label 
(1=1) ,  one  is  assigned  to  all  the  unspecified  bits  by  increasing  the 
index  I.   SHIFT  is  a  variable  which  represents  the  weight  assigned  to 
each  unspecified  bit  indicated  by  index  I. 

In  block  2,  the  value  of  USPFY  is  examined  and  if  it  is  larger 
than  M  (the  number  of  output  functions) ,  then  in  the  loop  which  consists 
of  blocks  3,  4,  5,  6  and  7,  the  values  of  output  functions  assigned  to 
vertex  B  (M  bits)  is  examined.   In  this  loop,  each  time  when  a  don't 
care  bit  is  encountered,  one  is  assigned  to  the  don't  care  bit  by  adding 
SHIFT  to  LBX.   Then,  in  the  loop  which  consists  of  blocks  8,  9  and  10, 
one  is  assigned  to  every  remaining  unspecified  bit.   If  the  value  of 
USPFY  does  not  exceed  M,  then  the  unspecified  bits  in  the  label  field 
of  vertex  B,  the  number  of  which  is  specified  in  input  parameter  USPFY, 
are  examined  and  one  is  assigned  to  every  unspecified  don't  care  bit. 

(9)   Subroutine  DCRMNT 

This  subroutine  has  five  parameters;  three  of  them  (NDCARE,  BWEIT 
and  M)  are  input  parameters  and  the  others  (INC,  LBXX)  are  input-output 
parameters.   In  this  subroutine,  parameter  LBXX  is  decreased  to  the 
next  possible  state.   As  shown  in  block  16  in  Fig.  3.3.12  (CMXL) ,  prior 
to  calling  this  subroutine,  LBXX  has  to  be  initialized  to  the  value 
stored  in  LBX. 
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Fig.  3.3.14   Flowchart  of  subroutine  ASIGN1. 
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As  can  be  seen  in  the  flowchart  (Fig.  3.3.15),  this  subroutine  is 
executed  in  the  same  way  as  subroutine  INCMNT.   The  only  exception  is 
blocks  7  and  8.   In  these  blocks,  BWEIT(J)  and  SHIFT  are  subtracted 
from  LBXX  instead  of  being  added. 

(10)   Subroutine  MPF 

Flowchart  is  shown  in  Fig.  3.3.16  and  the  following  is  the  defini- 
tion of  variables  used  in  this  subroutine. 

MXVTXL:   This  is  a  variable  which  stores  the  maximum  value, 

N 
(2  -1) ,  of  vertex  number  in  the  N-cube. 

NVTEXF :  This  is  a  variable  which  stores  the  number  of  vertices, 

(2      ) ,  in  the  Large-cube. 

RF— TT 
RSHIFT :   This  is  a  variable  which  stores  value  2 

LSHIFT :   This  is  a  variable  which  stores  value  2    . 

I_:   This  is  an  index  of  the  program  loop  in  Fig.  3.3.16  and  at  the 
same  time  I  represents  a  vertex  number  in  the  N-cube. 

.J:   This  is  a  variable  which  stores  the  vertex  number  of  a  vertex 
in  the  Large-cube  to  which  the  maximum  permissible  function  value  is  to 
be  assigned.   As  can  be  seen  in  Fig.  3.3.16,  when  vertex  I  in  the  N- 
cube  is  concerned,  the  corresponding  II-th  MPF  value  will  be  assigned 
to  vertex  J  in  the  Large-cube,  where  J  is  calculated  by  the  following 
formula  (see  Fig.  3.3.17); 

J  =  I  *  2II_1  +  LABEL  (1+1)  /  2RF~II+1 

RF— T  T 

=  I  *  LSHIFT  +  MNL  (1+1)  /  2  *  2 
-  I  *  LSHIFT  +  MNLX  /  2 


( 


START 
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IHC  *  INC  ♦  I 


LNOC  *   INC 


SHIFI  -  2 

J  ♦  1 


M-NDCARE 


LBXX  *  LBXX  -   BVEIT(J) 


RETURN 


LBXX  •■  LBXX   -   SHIFT 


I 


J  *  J  +   1 

SHIFT  »  2   *    SHIFT 


10 


INCX  •    INCX/2 


Fig.  3.3.15   Flowchart  of  subroutine  DCRMNT. 
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Fig.  3.3.16    Flowchart  of  subroutine  MTF. 
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Fig.  3.3.17   Input  vector  and  LABEL  of  vertex  I  in 
the  N-cube. 


In  block  1  (Fig.  3.3.16),  the  FUNC  field  of  the  Large-cube  is 
initialized  to  zero.   In  the  loop  which  consists  of  blocks  3,  4, . . . ,10, 
the  II-th  most  significant  bits  of  MNL  (1+1)  and  MXL  (1+1)  are  compared. 
If  the  resulting  MPF  value  is  "0,"  1  is  assigned  to  FUNC  (J+l) ,  if  the 
MPF  value  is  "1,"  2  is  assigned  to  FUNC  (J+l)  and  if  the  MPF  value  is 
"don't  care,"  0  is  assigned  to  FUNC  (J+l).   (In  the  last  case  nothing 
should  be  done  because  FUNC  is  already  initialized  to  zero  in  block  1.) 
We  should  notice  that  MNL  (1+1)  and  MXL  (1+1)  store  the  minimum  and  the 
maximum  labels  assigned  to  vertex  I  in  the  N-cube,  respectively,  and 
FUNC  (J+l)  stores  the  FUNC  field  of  vertex  J  in  the  Large-cube.  This 
loop  is  repeated  until  all  the  vertices  in  the  N-cube  are  exhausted. 


(11)   Subroutine  IMC 

As  described  in  Section  3.2,  this  subroutine  obtains  an  irredundant 
MOS  cell  configuration  from  the  maximum  permissible  function  in  the 
Large-cube.  This  subroutine  consists  of  six  steps  which  will  be 
described  below.   Before  describing  each  step  in  detail,  the  definitions 
of  variables  which  appear  in  these  steps  in  common  are  explained  first. 
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(The  definitions  of  variables  MXVTXL,  NVTEXF,  LSHIFT  and  RSHIFT  already 
appeared  in  the  description  of  subroutine  MPF.) 

INPD:  This  is  a  variable  which  stores  the  dimension  of  the  Large- 
cube  which  is  calculated  in  N  +  II  -  1.   This  variable  also  represents 

the  dimension  of  input  vectors  in  the  Large-cube. 

N 
NVTEXL:   This  is  a  variable  which  stores  the  number,  2  ,  of 

vertices  in  the  N-cube. 

MXVTXF :   This  is  a  variable  which  stores  the  maximum  value, 
2      -1,  of  vertex  number  in  the  Large-cube. 

ERROR:   This  is  a  logical  variable  which  stores  an  error  status 
detected  in  Step  6  of  subroutine  IMC. 

Step  1:   This  step  constructs  a  Large-cube  in  the  same  way  as  sub- 
routine NCUBE  constructs  a  N-cube.   The  flowchart  is  shown  in 
Fig.  3.3.18.   In  this  flowchart  variable  Y  corresponds  to  variable  X 
in  subroutine  NCUBE.   Other  variables  have  the  same  definitions  as  in 
subroutine  NCUBE. 

Step  2:   The  flowchart  for  step  2  is  shown  in  Fig.  3.3.19.   In  block  1, 
the  FUNC  field  of  the  vertex  with  weight  INPD  is  examined  and  if  it 
contains  don't  care,  "0"  is  assigned  to  the  vertex.   Then  all  the 
vertices  in  the  Large-cube  are  traversed  in  the  same  way  as  in  sub- 
routine CMNL  and  every  vertex  which  was  assigned  don't  care  in  subroutine 
MPF  is  assigned  "0"  or  "1"  in  the  following  way. 

As  in  subroutine  CMNL,  vertex  B  to  which  value  "0"  or  "1"  is  going 
to  be  assigned  is  pointed  by  PTRB.   The  values  assigned  to  all  the 
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Fig.  3.3.18   Flowchart  of  subroutine  IMC.   (Step  1) 
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Fig.  3.3.19   Flowchart  of  subroutine  IMC.   (Step  2) 
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vertices  which  are  connected  to  vertex  B  by  the  edges  from  these 
vertices  are  examined  in  the  loop  which  consists  of  blocks  7 , ... ,12. 
If  all  these  vertices  take  the  value  "0"  (1  or  3),  then  vertex  B  is 
assigned  "0"  by  storing  3  to  the  FUNC  field.   Otherwise  vertex  B  is 
assigned  "1"  by  storing  2  to  the  FUNC  field. 

In  order  to  get  an  irredundant  cover,  it  is  necessary  to  distin- 
guish the  original  "0"  which  has  already  been  assigned  in  subroutine 
MPF  from  the  "0"  which  is  assigned  in  this  step.  This  is  accomplished 
by  assigning  1  or  3  to  the  FUNC  field  of  the  Large-cube  according  to 
the  type  of  "0"  (1  for  original  "0"  and  "3"  for  the  "0"  which  is 
assigned  in  this  step).   In  the  following  flowcharts,  the  original  "0" 
is  expressed  as  "0*"  and  the  "0"  which  is  assigned  in  this  step  is 
expressed  as  "0." 

Step  3;  This  step  obtains  the  set  of  minimum  vectors.  The  flowchart 
for  step  3  is  shown  in  Fig.  3.3.20.   In  this  step,  NMINV  is  a  variable 
which  stores  the  number  of  minimum  vectors  and  MINV  is  an  array  which 
stores  the  minimum  vectors  in  binary  form. 

All  the  vertices  in  the  Large-cube  are  traversed  in  the  same  way 
as  in  subroutine  CMNL  and  every  vertex  to  which  zero  (both  "0"  and 
"0*"  are  included)  is  assigned  is  examined  to  determine  whether  the 
input  vector  assigned  to  the  vertex  belongs  to  the  set  of  minimum 
vectors  or  not  in  the  following  way. 

As  in  subroutine  CMNL,  a  vertex  B  under  examination  is  pointed  by 
PTRB.   The  values  assigned  to  all  the  vertices  which  are  connected  to 
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Fig.  3.3.20   Flowchart  of  subroutine  IMC.   (Step  3) 
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vertex  B  by  the  edges  from  vertex  B  to  these  vertices  are  examined  in 
the  loop  which  consists  of  blocks  6, ... ,11  and  if  all  these  vertices 
take  value  "1,"  the  input  vector  assigned  to  vertex  B  is  added  to  the 
set  of  minimum  vectors. 

After  the  implementation  of  step  3,  the  minimum  vectors  are 
obtained  in  array  MINV  and  the  number  of  the  minimum  vectors  are 
obtained  in  variable  NMINV. 

Step  4:   This  step  obtains  the  subset  of  the  set  of  minimum  vectors 
which  covers  all  the  vertices  with  "0*"  (original  zero)  and  this  sub- 
set is  called  a  semi-irredundant  subset.   Flowchart  is  shown  in 
Fig.  3.3.21. 

I_  is  an  index  which  points  to  a  minimum  vector  in  array  MINV. 

J  is  an  index  which  points  to  a  vertex  in  the  Large-cube. 
In  block  2,  the  first  minimum  vector  MINV(l)  is  picked  up  and  put 
into  the  semi-irredundant  subset.   In  the  loop  which  consists  of 
blocks  4,  5, . . . ,11,  the  FUNC  fields  of  all  the  vertices  in  the  Large- 
cube  are  examined  in  the  order  stored  in  the  array  FUNC.  The  input 
vector  of  each  vertex  with  "0*"  (original  zero)  is  compared  with 
MINV(l)  in  block  6  and  if  MINV(l)  covers  the  vertex  with  "0*,"  the 
LINK  field  of  the  vertex  is  increased  by  one.   (From  now  on,  since  we 
do  not  use  the  cube  structure  constructed  in  step  1,  LINK  field  is 
used  for  storing  the  degree  which  indicates  how  many  corresponding 
vertices  with  "0*"  are  covered  by  a  subset  of  the  set  of  minimum 
vectors.   The  LINK  field  is  initialized  to  zero  in  block  1.) 


1                STAJ1T 

1                       |l 

INITIALIZE 
LINK(NVTEXF) 

1                   V 

I  •  1 

3                     V 

REPEAT   •     .FALSE. 
J   •    1 

JX  -   J  -  1 
MINX  -  MINV(I) 
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Fig.  3.3.21   Flowchart  of  subroutine  IMC.   (Step  A) 
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After  exhausting  all  the  vertices  in  the  Large-cube,  if  there 
exists  a  vertex  with  "0*"  which  is  not  covered  at  all  (LINK  field  =  0), 
the  next  minimum  vector  is  picked  up  in  block  13  and  the  program  control 
returns  to  block  3.  The  same  procedure  is  repeated  until  all  the 
vertices  with  "0*"  are  covered  with  a  subset  of  the  set  of  minimum 
vectors. 

Since  the  set  of  minimum  vectors  is  constructed  such  that  it  covers 
all  the  zero's  in  Large-cube  (both  "0"  and  "0*"  are  included),  it  is 
obvious  that  the  semi-irredundant  subset  obtained  in  the  above  procedure 
can  cover  all  the  vertices  with  "0*." 

NSIRR  stores  the  number  of  input  vectors  in  the  semi-irredundant 
subset  obtained  in  the  above  procedure. 

Step  5;  This  step  obtains  an  irredundant  subset  from  the  semi-irredundant 
subset  obtained  in  step  A.   Flowchart  is  shown  in  Fig.  3.3.22.   In  this 
step,  I_  and  J_  are  defined  in  the  same  way  as  in  step  A. 

In  block  1,  the  first  minimum  vector  in  the  semi-irredundant  sub- 
set, MINV(l) ,  is  picked  up  and  is  tried  to  be  removed  from  the 
semi-irredundant  subset.  This  is  implemented  in  the  following  way. 
In  the  loop  which  consists  of  blocks  3,  A, ...  ,9,  the  FUNC  fields  of 
all  the  vertices  in  the  Large-cube  are  examined  and  the  input  vector 
of  each  vertex  with  "0*"  is  compared  with  MINV(l). 

If  MINV(l)  covers  the  vertex  with  M0*,M  the  LINK  field  of  the 
vertex  is  decreased  by  one  because  we  are  now  trying  to  remove  MINV(l) 
from  the  semi-irredundant  subset.   In  block  7,  the  value  of  the  result- 
ing LINK  field  is  examined  and  if  the  value  is  zero  (This  means  if  we 
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Fig.  3.3.22   Flowchart  of  subroutine  IMC.   (Step  5) 
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remove  MINV(l)  from  the  semi-irredundant  subset,  the  resulting  subset 
does  not  cover  all  the  vertices  with  "0*"  any  longer.) »  MINV(l)  cannot 
be  removed  from  the  semi-irredundant  subset.   The  LINK  fields  which 
has  been  decreased  in  the  loop  which  consists  of  blocks  3, ... ,9  have 
to  be  restored  in  the  loop  which  consists  of  blocks  11, ... ,16  by  in- 
creasing them  by  one.  After  the  restoring  operation,  the  next  minimum 
vector  is  picked  up  in  block  17 . 

After  exhausting  all  the  vertices  in  the  Large-cube  comparing  with 
MINV(l),  if  all  the  vertices  with  "0*"  are  still  covered  with  the 
semi-irredundant  subset  (There  exist  no  vertices  with  "0*"  whose  LINK 
fields  take  the  value  zero.),  MINV(l)  can  be  removed  from  the  semi- 
irredundant  subset  by  assigning  negative  value  to  MINV(l)  (block  10) . 
The  next  minimum  vector  is  picked  up  in  block  17  and  the  program  con- 
trol returns  to  block  2. 

The  above  process  is  repeated  until  all  the  input  vectors  except 
the  last  one  in  the  semi-irredundant  subset  are  exhausted.   The  result- 
ing subset  is  called  the  irredundant  subset  which  realizes  an  irredundant 
MOS  cell  configuration.   Finally,  in  block  19,  the  irredundant  MOS  cell 
configuration  is  printed  (The  output  format  is  described  in  detail  in 
Section  5.1.) . 

The  number  of  input  vectors  in  the  irredundant  subset  is  obtained 
in  a  variable  NIRR. 

Step  6;  This  step  stores  the  function  value  realized  by  the  MOS  cell 
obtained  in  step  5  to  the  II-th  most  significant  bit  of  the  LABEL  field 
in  the  N-cube.   Two  flowcharts  are  shown  in  Figs.  3.3.23  and  3.3.24. 
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Fig.  3.3.23    Flowchart  of  subroutine  IMC.   (Step  6A) 
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LABELCI+l)    ♦ 
LABEL(I+1)+RSHIFT 


Fig.  3.3.24   Flowchart  of  subroutine  IMC.   (Step  6B) 
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In  these  flowcharts,  I_  is  the  index  of  the  program  loop  and  at  the  same 
time  I_  represents  the  vertex  number  in  the  N-cube.  J[  is  the  variable 
which  stores  the  vertex  number  of  the  vertex  in  the  Large-cube  to  which 
the  maximum  permissible  function  (which  is  completely  specified)  is 
assigned. 

In  block  1  (Fig.  3.3.23),  the  dimension  of  the  irredundant  subset 
is  compared  with  the  dimension  of  the  set  of  minimum  vectors  and  if 
they  coincide,  the  function  realized  by  the  irredundant  MOS  cell  takes 
the  same  values  as  those  assigned  to  the  Large-cube  in  step  2.  There- 
fore, in  this  case,  the  function  value  which  has  been  obtained  at 
vertex  J  in  the  Large-cube  is  just  copied  into  the  II-th  most  signifi- 
cant bit  of  the  LABEL  field  of  vertex  I  in  the  N-cube.  As  described 
in  subroutine  MPF,  there  exists  the  following  relation  between  I  and  J; 

J  =  I  x  LSHIFT  +  LABEL  (1+1)  /  (2  *  RSHIFT) 

The  above  procedure  is  described  in  flowchart  in  Fig.  3.3.24. 

On  the  other  hand,  in  block  1  (Fig.  3.3.23),  if  they  take  dif- 
ferent values  (this  occurs  in  most  cases),  values  of  the  function 
realized  by  the  irredundant  MOS  cell  are  different  from  those  assigned 
to  the  Large-cube  in  Step  2.   Flowchart  for  this  case  is  shown  in 
Fig.  3.3.23.   In  this  flowchart,  after  determining  the  function  value 
realized  at  vertex  J  in  the  Large-cube  by  the  MOS  cell  obtained  in 
step  5  (this  is  accomplished  by  the  loop  which  consists  of  blocks  5,6, 
.  .  .  ,9) ,  it  is  stored  in  the  II-th  most  significant  bit  of  the  LABEL 
field  of  vertex  I  in  the  N-cube  (block  13) .   If  the  function  value  is 
zero  no  action  is  taken  because  zero  has  already  been  stored. 
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When  the  above  function  value  is  stored  in  a  bit  of  the  LABEL 
field  in  which  the  output  function  value  has  already  been  stored,  the 
two  function  values  are  compared.   If  they  do  not  coincide,  an  error 
flag  is  set  (blocks  14,15). 

After  storing  all  the  function  values  in  the  II-th  most  significant 
bit  of  LABEL  field,  the  status  of  an  error  flag  is  examined  in  block  12 
(Fig.  3.3.24)  and  if  an  error  flag  is  on,  an  error  message  and  the  con- 
tents of  array  LABEL  are  printed.  After  the  contents  of  array  LABEL 
is  printed  by  calling  subroutine  PRINT,  the  program  control  returns  to 
the  START  of  subroutine  MAIN  and  the  next  problem  is  read  in.   If  there 
is  no  error,  the  control  returns  to  subroutine  MAIN. 

(12)   Other  Subroutines 

(a)  Subroutine  COMPR  -  In  this  subroutine,  two  input  vectors  JX 
and  MINX  are  compared  bit  by  bit  and  if  every  bit  in  JX  is  larger  than 
or  equal  to  the  corresponding  bit  in  MINX,  the  program  control  takes 
the  normal  return.   Otherwise  the  control  is  transferred  to  the  state- 
ment number  shown  in  the  parameter  list. 

(b)  Subroutines  ASFUN1  and  ASFUN2  -  These  subroutines  are 
explained  in  Section  3.2.   Flowcharts  are  shown  in  Fig.  3.3.25  and 
Fig.  3.3.26. 
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Fig.  3.3.25   Flowchart  of  subroutine  ASFUN1 


START 


69 


YES 


INITIALIZE 
KUNC(NVTEXF) 


1  •  0 


MNLX   •    MNI.(I+1)/RSHIFT 
J    ■     1    *   1.SHIFT   +  MMLX/2 


YES 


FUNC(J+1)    •>    1    ("0") 


I   -   I   +   1 


NO 


FUNC(.J+1)    -    2    ("1") 


RETURN 


Fig.  3.3.26   Flowchart  of  subroutine  ASFUN2, 
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4.  INPUT  DATA  SETUP 

4.1  Input  Data  Card  Format 

For  each  separate  problem,  a  set  of  input  data  cards  must  be  sub- 
mitted which  consist  of  the  following. 

(i)   <  parameter  card  > 

(ii)  <  output  function  card  >  s 
(i)  will  always  consist  of  only  a  single  card,  but  (ii)  may  consist  of 
more  than  one  card.   The  format  of  each  data  card  is  explained  in  the 
following. 

(i)   <  parameter  card  > 

This  card  specifies  the  number  of  external  variables,  N,  and  the 
number  of  output  functions,  M,  for  a  given  problem.   The  parameter  N 
and  M  are  specified  in  the  following  two  fields. 

Cols.  1-2:   This  field  specifies  an  integer,  N,  which  is  right  justified 
Cols.  3-5:   This  field  specifies  an  integer,  M,  which  is  right  justified, 

(ii)   <  output  function  card  >  s 

Although  output  functions  are  submitted  to  this  program  in  the 
truth  table  form,  the  input  vector  for  each  output  function  value  is 
implicitly  specified  by  the  order  in  which  the  function  values  are 
specified  in  each  output  function  card(s).   Depending  on  the  function 
value  (one,  zero  or  don't  care),  a  different  character  ("1,"  "0"  or  "*") 
is  punched  on  the  corresponding  column. 

The  truth  table  with  N  external  variables  and  M  output  functions  is 
shown  in  Fig.  4.1.   This  truth  table  is  submitted  to  the  program  DIMN 
in  the  fo] lowing  way. 


71 


xl  x2   •  •  •   xn 


LM 


0   0   .  .  .   0 
0   0   .  .  .   1 


1   1 


1   1 


.2N-2 


.2N-1 


;2N_2 
2 

.2N-1 


"M 


M 


:2N-2 

"M 

:2N-1 
'M 


Fig.  4.1   Truth  table  with  N  external  variables  and 
M  output  functions. 


Each  output  function  is  specified  on  one  or  several  output  function 
cards  depending  upon  the  number  of  external  variables,  N.   Since  each 
function  value  is  specified  in  one  column  in  the  output  function  card(s), 

if  the  number  of  external  variables,  N,  is  larger  than  or  equal  to  7, 

N 
the  number  of  input  vectors,  2  ,  will  exceed  the  number  of  columns  in 

one  card  and  consequently  two  or  more  output  function  cards  are  required 

in  order  to  specify  one  output  function.   In  fact,  the  number  of  cards 

f~2N  "I 

needed  to  specify  one  output  function  is  equal  to   —   .   (That  is, 

oU 

^  9 

if  N=9,  then  2  =  512.   This  means  that  we  need  7  cards  to  specify  one 

output  function.) 
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In  Fig.  4.1,  f,  ,  the  first  function  value  of  function  f   is  speci- 
fied in  the  first  column  of  the  output  function  card(s)  which  specify 
function  f  .   f .  ,  the  second  function  value  of  function  f.  is  specified 
in  the  second  column  and  so  forth.   After  specifying  all  the  function 
values  of  function  f1 ,  function  f ^  is  specified  on  separate  output 
function  card(s).   The  above  process  is  repeated  until  M  output  func- 
tions are  specified  on  M  separate  sets  of  output  function  card(s). 
Since  blank  column  terminates  one  output  function  specification,  (the 
program  DIMN  interprets  the  blank  character  as  the  termination  sign  of 
one  output  function  specification)  the  blank  character  should  not  be 
inserted  among  function  specification  except  at  the  end  of  each  func- 
tion specification. 

In  Fig.  4.2,  input  data  cards  for  one  problem  is  shown  and  in 


S 


input 

data 

cards 


c 


output 

function 

cards 


c 


7 


c 


< 


7 


parameter  card  > 


output  function 
specification  for  f. 


output  function 
specification  for  f~ 


output  function 
specification  for  f 


M 


Fig.  A. 2   The  input  data  cards  for  one  problem. 
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Fig.  4.3,  an  input  card  sequence  for  the  execution  of  a  typical  DIMN 
problem  is  shown.   We  can  put  as  many  problems  as  we  want  in  the  fields 
of  these  input  data  cards.   If  there  exists  a  format  error  in  input 
data,  the  following  error  message  is  printed  out  and  the  program  execu- 
tion halts. 

"INPUT  ERROR  IN  DATA  CARD  I  =  ***  J  =  ***  K  -  ***" 

This  message  means  that  an  error  occurred  at  the  K-th  column  of  the 
J-th  output  function  card  which  specifies  the  I-th  output  function. 

4.2  Restriction  on  Problem  Size 

In  order  to  pack  problems  into  a  finite  amount  of  space,  some  re- 
strictions on  the  size  of  an  acceptable  problem  are  required.   In  any 
case,  the  sum  of  the  number  of  external  variables,  N  and  the  number  of 
output  functions,  M  must  not  exceed  11  and  the  value  of  M  must  not 
exceed  4.  The  relation  between  the  maximum  number  of  external  vari- 
ables and  the  maximum  number  of  output  functions  which  the  program 
DIMN  can  handle  is  shown  below. 

Max.  no.  of  external  variables      Max.  no.  of  outputs 

10  1 

9  2 

8  3 

7  4 

If  the  above  restriction  is  violated,  the  following  error  message 
is  printed  and  the  program  execution  halts. 

"INPUT  ERROR  IN  PARAMETER  CARD  N  =  ***  M  =  ***" 
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/*  ID   <  ID  card  information  > 

/*  ID  REGION  =  200K,  TIME  =  (00.30),  LINE  =  03000 

//   EXEC  FORTLDGO,  REGION.  GO  -  200K 


FORTRAN  source  program 


/* 

//   GO.  SYSIN  DD  * 


input 

data 

cards 


< 


C 


L 


f 


c 


L 


<  parameter  card  > 


<  output-function-card  >  s 


parameter  card 


<  parameter  card  > 


the  first  problem 


the  second  problem 


the  last  problem 


/* 


Fig,  k  . 3   Input  card  sequence  for  the  execution  of  a  typical 
DIMN  problem. 


75 


In  this  error  message,  N  is  the  number  of  external  variables  and  M  is 
the  number  of  output  functions  on  a  parameter  card.   These  limitations 
are  essentially  imposed  by  the  array  sizes  in  the  programs  presently 
written.  To  loosen  the  restrictions  is  mainly  a  task  of  increasing 
array  dimensions  appropriately. 

4.3  Example  of  Input  Data  Setup 

The  following  examples  show  how  the  input  data  is  set  up  for  the 
typical  DIMN  problems. 

Example  1:  The  input  data  setup  for  the  network  with  three  external 
variables  (N=3)  and  one  completely  specified  output  function  (M=l)  is 
shown  in  Fig.  4.5.  The  truth  table  for  the  function  is  shown  in 
Fig.  4.4. 

Example  2:  The  input  data  setup  for  the  network  with  three  external 
variables  (N=3)  and  two  incompletely  specified  output  functions  (M=2) 
is  shown  in  Fig.  4.7.  The  truth  table  for  the  functions  is  shown  in 
Fig.  4.6. 
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X1X2 

X3 

fl 

0     0 

0 

0 

0     0 

1 

0 

0     1 

0 

1 

0     1 

1 

0 

1     0 

0 

0 

1     0 

1 

0 

1    1 

0 

1 

1    1 

1 

1 

Fig.  4.4   Truth  table  for  Example  1. 


column  no. 


000000000111 

123456789012 


output  function 
card  (f  ) 


parameter  card 


Fig.  4.5   Possible  setup  of  data  cards  to  specify  the 
problem  given  in  Example  1. 
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x_   x_   x« 


£1 

f2 

1 

* 

* 

0 

0 

1 

1 

* 

* 

0 

* 

* 

0 

1 

1 

1 

Fig.  4.6   Truth  table  for  Example  2. 


column  no. 


000000000111 

123456789012 


output  function 
card  for  f„ 


output  function 
card  for  f. 


parameter  card 


Fig.  4.7   Possible  setup  of  data  cards  to  specify  the 
problem  in  Example  2 . 
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5.  OUTPUT  OF  PROCEDURE  DIMN 

This  chapter  describes  the  output  of  program  DIMN.   In  Section  5.1, 
the  output  format  is  shown  for  a  typical  DIMN  problem  and  in  Section  5.2, 
MOS  networks  obtained  by  program  DIMN  are  compared  with  the  correspond- 
ing MOS  networks  obtained  by  previously  developed  algorithm  for  several 
4  variable  functions. 


5.1  Output  Format 

Fig.  5.1  shows  the  printout  obtained  by  program  DIMN  for  a  typical 
example.   In  this  printout,  the  subscripted  variables  XI,  X2  -  -  -  - 
represent  the  external  input  variables  connected  to  driver  MOS  FETs  and 
the  subscripted  variables  Ul,  U2,  -  -  -  -  represent  the  outputs  of  MOS 
cells  connected  to  the  inputs  of  other  driver  MOS  FETs. 

This  printout  shows  the  number  of  external  variables  and  the 
number  of  output  functions;   two  and  three,  respectively.   Then  the 
two  output  functions  are  printed  in  the  truth  table  form  where  input 
vectors  are  implicitly  specified  by  the  order  in  which  the  function 
values  appear. 

In  the  printout  of  MOS  cell  configurations,  only  driver  MOS  cells 
are  printed.   Fig.  5.2  shows  the  network  obtained  from  the  MOS  cell 
configuration  printed  in  Fig.  5.1.   In  this  network,  factoring  of 
literals  in  switching  expressions  is  done  by  hand. 

Following  the  above  MOS  cell  configurations,  some  statistics  which 
are  derived  from  the  obtained  network  are  printed.   In  this  example, 
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***************************************** 

*  * 

*  DESIGN  OF  IRREDUNDANT  MOS  NETWORK    * 

*  * 
***************************************** 

X  a  EXTERNAL  VARIABLE 
U  =  OUTPUT  OF  MOS  CELL 

***************************************** 


NUMBER  OF  EXTERNAL  VARIABLES  = 
NUMBER  OF  OUTPUT  FUNCTIONS 

FUNCTION   1 
10010101 

FUNCTION   2 
00100111 

NETWORK  CONFIGURATION 

MOS  CELL  I   0 X2 — X3 0 


MOS  CELL  2 


I — X2 — LI— 
0 —  !  —  0 

I  —  XI — X3 — 


MOS  CELL  3   0 — 


— 0 


MOS  CELL  4   0 U2 0 


NUM8ER  OF  MOS  CELLS  =  4 

NUMBER  OF  MOS  FETS   =  14 
(  WITHOUT  FACTORING  ) 

ELAPSEO  TIME  =    0.12  SEC 


Fig.  5.1   Printout  obtained  by  program  DIMN. 
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the  number  of  MOS  cells  is  four,  the  number  of  driver  MOS  FETs  is  14 
(without  factoring)  and  elapsed  time  for  constructing  the  network  is 
0.14  sec.   As  a  timing  subroutine,  STEPZ  in  FORTRAN  system  subroutines 
is  used. 

5.2  Networks  Obtained  by  Program  DIMN 

The  networks  obtained  by  program  DIMN  are  shown  in  Appendix  A  for 
several  three  and  four  variable  functions.   Each  page  in  Appendix  A 
contains  the  output  of  program  DIMN  for  one  problem  in  upper  half  and 
the  resulting  network  obtained  from  the  above  computer  output  in  lower 
half.   In  the  resulting  networks,  literals  in  switching  expressions  are 
factored  by  hand. 

Comparing  the  four  variable  networks  obtained  in  Appendix  A  with 
the  corresponding  four  variable  networks  obtained  in  [6]  based  on 
previously  developed  Liu's  algorithm,  we  observe  that  significant 
improvement  with  respect  to  the  number  of  driver  MOS  FETs  can  be 
achieved  by  Lai's  algorithm  on  which  program  DIMN  is  based. 

The  number  of  MOS  FETs  obtained  by  two  different  algorithms  for 
twelve  4-variable  single-output  functions  are  compared  in  Table  5.1. 

In  this  table,  functions  are  represented  in  hexadecimal  form. 
For  example,  function  f  =  000000000000  1011  is  represented  in  hexa- 
decimal form  f  =  000B.  This  table  shows  that  the  average  number  of 
MOS  FETs  is  improved  by  a  factor  of  1.36  (from  12.9  to  9.5)  by 
Algorithm  DIMN.   In  the  following  we  can  show  one  typical  example  for 
the  function  f,  =  000B  (hexadecimal  representation)  with  4  external 
variables. 
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Function  in 
hexadecimal 

Number  of  driver  MOS  FETs  in  Network 

obtained  by  Liu's  Algo. 

i 

obtained  by  Lai's  Algp. 

0001 

5 

5 

0002 

6 

5 

0006 

9 

7 

0007 

5 

5 

0008 

7 

5 

0009 

13 

8 

000B 

12 

6 

0016 

12 

10 

0069 

21 

12 

0117 

9 

9 

6996 

28 

21 

9669 

28 

21 

Total  FETs 

155 

114 

Average  FETs 
in  one  Network. 

12.9 

9.5 

Table  5.1   Number  of  driver  MOS  FETs  in  networks  obtained  by 
two  different  algorithms. 


Program  DIMN  is  based  on  Lai's  algorithm. 
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The  network  obtained  by  program  DIMN  and  the  network  obtained  based 
on  Liu's  algorithm  are  shown  in  Fig.  5.3  and  Fig.  5.4,  respectively.   As 
can  be  seen  in  these  figures,  the  number  of  driver  MOS  FETs  for  this 
particular  example  can  be  reduced  to  one  half  by  applying  Algorithm  DIMN, 
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Fig.  5.3   Network  obtained  by  program  DIMN  for  f..  =  000B; 
3  load  MOS  FETs  and  6  drivers. 


HC 


x2  <H 


O  f . 


Fig.  5.4   Network  obtained  applying  previously  developed 
algorithm  for  f   =  000B;  3  load  MOS  FETs  and 
12  drivers. 
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6.  CONCLUSION 

As  discussed  in  the  introduction,  Lai's  algorithm  extended  the 
previously  developed  algorithms  and  made  it  possible  to  obtain  irre- 
dundant  MOS  networks  with  minimum  number  of  negative  gates  (MOS 
cells) .   Although  Algorithm  DIMN  does  not  guarantee  exhaustion  of  all 
possible  irredundant  MOS  networks  for  an  arbitrary  given  function, 
Algorithm  DIMN  is  so  far  the  only  existing  procedure  which  can  obtain 
irredundant  MOS  networks  efficiently.   Our  computational  experience 
with  Algorithm  DIMN  implemented  as  program  DIMN  shows  the  efficiency 
of  this  algorithm. 
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APPENDIX  A 
Networks  Obtained  by  Program  DIMN 


***************************************** 

*  * 

*  DESIGN  OF  IRREOUNDANT  MOS  NETWORK    * 

*  * 
***************************************** 

X  =  EXTERNAL  VARIABLE 
U  *  OUTPUT  OF  MOS  CELL 

***************************************** 
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NUMBER  OF  EXTFRNAL  VARIABLES 
NUMBER  OF  OUTPUT  FUNCTIONS 

FUNCTION   1 
OOlOOOll 

NETWORK  CONFIGURATION 

MOS  CELL  1   0- X3 0 


MOS  CELL  2   0 — 


I  —  X2 — Ui — 
— XI — X2 — 


—  0 


MOS  CELL  3   0- U2 0 


NUMBER  OF  MOS  CELLS  =  3 
NUMBER  OF  MOS  FETS  =  6 
{  WITHOUT  FACTORING  ) 


ELAPSED  TIME  - 


0.07  SEC 


x3  ° — If 


***************************************** 

*  * 

*  DESIGN  OF  IRREDUNDANT  MOS  NETWORK    * 

***************************************** 

X  =  EXTERNAL  VARIABLE 
U  =  OUTPUT  OF  MOS  CELL 

***************************************** 
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NUMBER    OF    EXTERNAL    VARIABLES 
NUMBER    OF    OUTPUT    FUNCTIONS 

FUNCTION       1 

0*10***1 

NETWORK  CONFIGURATION 

MOS  CELL  1   0 X3 0 


MOS  CELL  2   0 — 


—  X2 — Ul- 
— XI 


I  —  0 


MOS  CELL  3   0 U2 0 

NUMBER  OF  MOS  CELLS  =  3 
NUMBER  OF  MOS  FETS  =  5 
(  WITHOUT  FACTORING  ) 

ELAPSED  TIME  =    0.07  SEC 


H 


4 


<—\ 


JH"i 


X 


Mr 


ii O  f 


c 


♦  * 

*  DESIGN  OF  IRREDUNDANT  MOS  NETWORK    * 

***************************************** 

X  =  EXTERNAL  VARIABLE 
U  =  OUTPUT  OF  MOS  CELL 
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NUMBER  OF  EXTERNAL  VARIABLES 
NUMBER  OF  OUTPUT  FUNCTIONS 

FUNCTION   1 
1*01**01 

FUNCTION   2 
*01*0*11 


NETWORK  CONFIGURATION 

MOS  CELL  1   0 X2 — X3 0 

MOS  CELL  2   0 X2 — Ul 0 

MOS  CELL  3   0 Ul — U2 0 


NUMBER  OF  MOS  CELLS  =  3 
NUMBER  OF  MOS  FETS  =  6 
(  WITHOUT  FACTORING  ) 


ELAPSED  TIME  = 


0.05  SEC 


x20— | 
x_0 1 
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***************************************** 

*  * 

*  DESIGN    OF    IRREDUNDANT    MOS    NETWORK         * 

*  * 
******  *********************************** 

X  =  EXTERNAL  VARIABLE 
U  =  OUTPUT  OF  MOS  CELL 

***************************************** 


NUMBER  OF  EXTERNAL  VARIABLES  =   4 
NUMBER  OF  OUTPUT  FUNCTIONS    =   1 

FUNCTION   I 
OOOOOOOOOOOOOOOl 

NETWORK  CONFIGURATION 

MOS  CELL  1   0 XI — X2 — XB — X4 0 


MOS  CELL  2   0- 


— Ul 0 


NUMBER  OF  MOS  CELLS  =  2 
NUMBER  OF  MOS  FETS  *  5 
(  WITHOUT  FACTORING  ) 


ELAPSED  TIME  = 


0.03  SEC 


xi  o—| 
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***************************************** 

*  * 

*  DESIGN  OF  IRREDUNDANT  MOS  NETWORK    * 

*  * 
***************************************** 

X  =  EXTERNAL  VARIABLE 
U  =  OUTPUT  OF  MOS  CELL 

***************************************** 


NUMBER    OF    EXTERNAL    VARIABLES    =       4 
NUMBER    OF    OUTPUT    FUNCTIONS         =       I 

FUNCTION       1 
OOOOOOOOOOOOOOIO 

NETWORK    CONFIGURATION 

MOS    CELL    1       0 XL — X2 — X3 0 


MOS    CELL    2       0 — 


— Ul — 
—  XA — 


—  0 


NUMBER  OF  MOS  CELLS  *  2 
NUMBER  OF  MOS  FFTS  =  5 
(    WITHOUT    FACTORING    ) 

ELAPSED    TIME    =         0.06  SEC 


O  f 


***************************************** 

*  * 

*  DESIGN  OF  IRREDUNDANT  MOS  NETWORK    * 

*  * 
**************** ******* *******  *********** 
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X  =  EXTERNAL  VARIABLE 
U  =  OUTPUT  OF  MOS  CELL 

***************************************** 


NUMBER  OF  EXTERNAL  VARIABLES  =   4 
NUMBER  OF  OUTPUT  FUNCTIONS    =   i 


FUNCTION   1 
0000000000000110 


NETWORK  CONFIGURATION 


MOS  CELL  1   0 — I 


—  Xi — X2 — X4 — 
— XI — X2 — X3 — 


I — 0 


MOS  CELL  2   0 — 


—  X3 — X4- 


—  0 


NUMBER  OF 
NUMBER  OF 
(  WITHOUT 


MOS  CELLS  =»  2 
MOS  FETS  =  9 
FACTORING  ) 


ELAPSED  TIME  -    0.06  SEC 


***************************************** 

*  * 

*  DESIGN  OF  IRREOUNDANT  MOS  NETWORK    * 

*  * 
***************************************** 

X  =  EXTERNAL  VARIABLE 
U  =  OUTPUT  OF  MOS  CELL 

***************************************** 
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NUMBER  OF  EXTERNAL  VARIABLES 
NUMBER  OF  OUTPUT  FUNCTIONS 

FUNCTION   I 
OOOOOOOOOOOOOlli 

NETWORK  CONFIGURATION 


4 
1 


MOS  CELL  1   0 — 


I  —  XI — X2 — X4 — 
'  —  XI — X2 — X3 — 


— 0 


MOS  CELL  2   0 Ul 0 


NUMBER  OF  MOS  CELLS  =  2 
NUMBER  OF  MOS  FETS  ■  7 
(  WITHOUT  FACTORING  ) 

ELAPSED  TIME  -    0.07  SEC 


***************************************** 

*  * 

*  DESIGN  OF  IRREDUNDANT  MOS  NETWORK    * 

*  * 
***************************************** 

X  =  EXTERNAL  VARIABLE 
U  =  OUTPUT  OF  HOS  CELL 

***************************************** 


NUMBER  OF  EXTERNAL  VARIABLES  =   4 
NUMBER  OF  OUTPUT  FUNCTIONS    =   1 
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FUNCTION       1 
OOOOOOOOOOOOIOOO 


NETWORK   CONFIGURATION 


MOS    CELL    1 


XI — X2 0 


MOS    CELL    2       0 — 


•Ul — 
X4 — 
X3 — 


—  0 


NUMBER    OF    MOS    CELLS    ■  2 

NUMBER    OF    MOS    FETS       =  5 
I    WITHOUT    FACTORING    ) 

ELAPSED    TIME    =         0.09  SEC 


xL0 1 


-Of 


sjHEH 


*    DESIGN  OF  IRREDUNDANT  MQS  NETWORK    * 

X  =  EXTERNAL  VARIABLE 
U  =  OUTPUT  OF  MOS  CELL 

***************************************** 
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NUMBER  OF  EXTERNAL  VARIABLES  = 
NUMBER  OF  OUTPUT  FUNCTIONS 

FUNCTION   1 
OOOOOOOOOOOOlOOl 

NETWORK  CONFIGURATION 

MOS  CELL  1   0 XI — X2 0 

MOS  CELL  2   0 X3 — X4 0 


4 
1 


MOS  CELL  3   0 — 


—  XA — U2- 

—  X3 — U2- 
— Ul 


— 0 


NUMBER  OF  MOS  CELLS  =  3 

NUMBER  OF  MOS  FETS   =  9 

(  WITHOUT  FACTCRING  I  . 

ELAPSED  TIME  =    0.15  SEC 


Of 


***************************************** 

*  * 

*  DESIGN  OF  IRREDUNDANT  MOS  NETWORK    * 

*  * 
***************************************** 

X  =  EXTERNAL  VARIABLE 
U  =  OUTPUT  OF  MOS  CELL 

***************************************** 
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NUMBER  OF  EXTERNAL  VARIABLES 
NUMBER  OF  OUTPUT  FUNCTIONS 

FUNCTION   1 
OOOOOOOOOOOOiOll 


NETWORK  CONFIGURATION 


MOS  CELL  I   0 XI — X2- 


MOS  CELL  2   0 X3 0 


MOS  CELL  3   0 — 


— X4 — U2 — 


—  0 


NUMBER  OF  MOS  CELLS  =  3 
NUMBER  OF  MOS  FETS  =  6 
(  WITHOUT  FACTORING  ) 

ELAPSED  TIME  =    0.15  SEC 


5  m 


H 


t 


_.    x4  ° II 


b 


■Of 


H 


T 


*  DESIGN  OF  IRREDUNDANT  MOS  NETWORK    * 

*  * 
***************************************** 

X  =  EXTERNAL  VARIABLE 
U  =  OUTPUT  OF  MOS  CELL 

***************************************** 
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NUMBER    OF    EXTERNAL    VARIABLES 
NUMBER    OF    OUTPUT    FUNCTICNS 

FUNCTION       1 
0000000000010110 


A 

1 


NETWORK  CONFIGURATION 


MOS  CELL  1   0 — 


— XI — X3 — X4- 

—  XI — X2 — X4- 

—  Xi — X2 — X3- 


—  0 


MOS  CELL  2   0 — 


—  X2 — X3 — X4~ 


—  0 


NUMBER  OF  MOS  CELLS  =  2 
NUMBER  OF  MOS  FETS  =  13 
(  WITHOUT  FACTORING  ) 


ELAPSED  TIME  = 


0.09  SEC 


O  f 


***************************************** 

♦  * 

*  DESIGN  OF  IRREDUNOANT  MOS  NETWORK    * 

***************************************** 

X  ■  EXTERNAL  VARIABLE 
U  =  OUTPUT  OF  MOS  CELL 

***************************************** 
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NUMBER  OF  EXTERNAL  VARIABLES  * 
NUMBER  OF  OUTPUT  FUNCTIONS 

FUNCTION   I 
00OOOOOOO11O1O01 

NETWORK  CONFIGURATION 

MOS  CELL  1   0 X3 — X4 0 


4 
1 


MOS  CELL  2 


XI — X« — Ul — 
XI — X3 — Ul  — 
XI — X2 


— 0 


MOS  CELL  3   0 — 


—  X2 — X4 — Ul 

—  X2 — X3 — Ul 


— 0 


NUMBER  OF  MOS  CELLS  *   3 

NUMBER  OF  MOS  FETS   *  17 
(  WITHOUT  FACTORING  ) 

ELAPSEO  TIME  =    0.19  SEC 


***************************************** 

*  * 

*  DESIGN  OF  IRREDUNDANT  MOS  NETWORK    * 

*  * 
***************************************** 

X  *  EXTERNAL  VARIABLE 
U  =  OUTPUT  OF  MOS  CELL 

***************************************** 
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NUMBER  OF  EXTERNAL  VARIABLES 
NUMBER  OF  OUTPUT  FUNCTIONS 

FUNCTION   1 
OOOOOOOlOOOlOlil 

NETWORK  CONFIGURATION 


4 
1 


MOS  CELL  1   0— 


—  X2 — X3 — XV 

—  Xi  —  X3 — X4- 

—  XI — X2 — XV 

—  XI  —  X2 — X3- 


— 0 


MOS  CELL  2   0 Ul 


NUMBER  OF  MOS  CELLS  =  2 
NUMBER  OF  MOS  FETS  =  13 
{  WITHOUT  FACTORING  ) 


ELAPSED  TIME  = 


0.07    SEC 


H 


5 


*3°— IL     xi°-l 
xa°HL     X2H 


O  f 


***************************************** 

*  * 

*  DESIGN  OF  IRREDUNDANT  MOS  NETWORK    * 

*  * 
***************************************** 

X  =  EXTERNAL  VARIABLE 
U  =  OUTPUT  OF  MOS  CELL 

***************************************** 
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NUMBER  OF  EXTERNAL  VARIABLES 
NUMBER  OF  OUTPUT  FUNCTIONS 

FUNCTION   I 
O11O1O011OO1011O 


4 

I 


NETWORK  CONFIGURATION 


MOS  CELL  1   0 — 


—  X4 — 

—  X3— 


—  0 


MOS  CELL  2   0— 


X2 — X3 — X4 — 
XI — X3 — X4 — 

X2 — Ul 

XI  — Ul 

Xi — X2 


—  0 


MOS  CELL  3   0 — 


—  XI — X2 — X3— X4- 

—  X3 — X4 — U2 

—  XI — X2 — Ul 

~  U I — U  2 

— X2 — U2 

—  XI — U2 


—  0 


NUMBER  OF  MOS  CELLS  =  3 

NUMBER  OF  MOS  FETS   =  30 
I  WITHOUT  FACTORING  ) 

ELAPSED  TIME  =    0,25  SEC 


H 


4 


£ 


\ 


><H 


***************************************** 

*  * 

*  DESIGN  OF  IRREOUNDANT  MOS  NETWORK    * 

**** ************************************* 

X  =  EXTERNAL  VARIABLE 
U  =  OUTPUT  OF  MOS  CELL 

***************************************** 


1 


NUMBER    OF    EXTERNAL    VARIABLES    = 
NUMBER    OF    OUTPUT    FUNCTIONS 

FUNCTION       1 
lOOiOllOOUOlOOl 

NETWORK    CONFIGURATION 

MOS    CELL     i       0 X3—  X4 0 

■X2 — X4 — Ul — 
X2 — X3 — Ul — 
■XI — X4 — Ul  — 
-XI — X3 — Ul — 
-XI — X2 


—  XI — X2 — XV—  Ul- 

—  XI — X2 — X3  —  Ul- 

—  X4— Ul — U2 

—  X3— Ul — U2 

— X2 — U2 

— XI — U2 


MOS  CELL  2   0 — 


—  0 


MOS  CELL  3   0 — 


—  0 


NUMBER  OF  MOS  CELLS  =   3 

NUMBER  OF  MOS  FETS   =  34 
(  WITHOUT  FACTORING  I 

ELAPSED  TIME  =    0.24  SEC 
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i  *i^i  x2°irii 


*HM 


H 


5 


i 


foffi  1  ^ 1  >^ 


103 


APPENDIX  B 
Program  Listing 


10/ 


C  THIS  PROGRAM  IS  THE  INPLEMENTATION  OF  ALGORITHM  DIMN  <  DESIGN  OF 

f  IRREDUNDANT  MOS  NETWORK  ). 

C  THIS  PROGRAM  IS  APPLICABLE  TO  THE  NETWORKS  WITH  MULTIPLE 

C  '   INCOMPLETELY  SPECIFIED  OUTPUTS. 

C  THE  RFLATICN  BETWEEN  THE  MAXIMUM  NUMBER  OF  EXTERNAL  VARIABLES  AND 

f.  THE  MAXIMUM  NUMBER  OF  OUTPUTS  WHICH  THIS  PROGRAM  CAN  HANDLE  IS 

C  SHOWEN  BELOW. 

C 

C  MAX.  EXTERNAL  VARIABLES        MAX.  OUTPUTS 

C 

C  10  1 

r  -9  ,2 

C  8  3 

C  7  4 

C 

C    THIS  PROGRAM  IS  CONSTRUCTED  WITH  THE  FOLLOWING  MAIN  PUBROUTINES  AND 

C    THE  OTHER  SMALL  SUBROUTINES. 

C 

C      SUBROUTINE  NAME     OPERATION 

C 

C  NCUBE  CONSTRUCT    N-CUBE    ACORDING    TO    THE    NUMBER    OF 

C  EXTERNAL    VARIABLES. 

C  INPUT  READ    DATA    INTO   N-CUBE    (     LABELt    DCARE    ). 

C  CMNL  IMPLEMENT    CONDITIONAL    MINIMUM    LABELING. 

C  CMXL     '  IMPLEMENT    CONDITIONAL    MAXIMUM    LABELING. 

C  MPF  ORTAIN    MAXIMUM    PERMISSIBLE    FUNCTION    FROM 

C  MINIMUM    LABEL    AND    MAXIMUM    LABEL. 

C  IMC  OBTAIN    IRREDUNDANT    MOS    CELL    CONFIGURATION 

C  FROM    MAXIMUM    PERMISSIBLE    FUNCTION. 

IMPLICIT  INTEGER!  A-Z  ) 

INTEGER+2   LABEL   tDCARE   ,MNL     ,MXL     , CHAIN   ,FUNC    ,LINK 
1  ,MINV 

COMMON       II       ,N        ,M        ,RF 

1  ,!ABEL(1024)     ,DCARE(1024)     ,MNL(1024)        ,MXL(1024I 

2  ,CHAIN(1024)     tSTARTLfll)      ,ENDL(11)         ,FUNC(16384) 

3  ,LINK( 16384)      ,MINV(  16384)      ,STARTF(15)       ,ENDF<15) 
LOGICAL   CS 

REAL  UTIME 
C    II  POINTER  WHICH  INDICATES  THE  CURRENT  STEP. 

C    N  NUMBER  OF  EXTERNAL  VARIABLES. 

C    M  NUMBER  OF  OUTPUT  FUNCTIONS. 

C    PF  MINIMUM  NUMBER  OF  MOS  CELLS 

C    CS  IS  SET  IN  SUBROUTINE-INPUT  IF  THE  GIVEN  OUTPUT 

C  FUNCTIONS  ARE  COMPLETELY  SPECIFIED. 

C    UTI"E  STORES  THE  TIME  ELAPSED  FOR  CONSTRUCTING  NETWORK. 

C    TFFT  STORES  TOTAL  NUMBER  OF  DRIVER  MOS  FETS  IN  NETWORK. 

C    LABEL(1024)      STORES  LABEL  AND  FUNCTION  VALUE  WHICH  IS  ASSIGNED 
r.  TO  EACH  VERTFX  IN  N-CUBE. 

C    OCARE(1024)       STORES  THE  DONTCARE  BIT  OF  OUTPUT  FUNCTIONS. 
C    MNL   (1024)       STORES  THE  LABEL  ASSIGNED  TO  EACH  VERTEX  BY 
C  CONDITIONAL  MINIMUM  LABELING. 

MXL   (1024)      STORES  THE  LABEL  ASSIGNED  TO  EACH  VERTEX  BY 
C  CONDITIONAL  MAXIMUM  LABELING. 

C    CHAINU024)       STOKES  THE  LINK  TO  THE  NEXT  VERTEX  WITH  THE  SAME 
C  WEIGHT  IN  N-CUBE. 

STARTL(ll)        POINTER  TO  THE  FIRST  VERTEX  WITH  THE  SAME  WEIGHT 
C  IN  N-CUBE. 

C    ENDL   (11)        POINTFR  TO  THE  LAST  VERTEX  WITH  THE  SAME  WEIGHT  IN 
C  N-CURE. 

FUNCl 163841       STURES  MAXIMUM  PERMISSIBLE  FUNCTION  ASSIGNED  TO  EACH 
C  VERTEX  IN  LAPGE  CUBE  FOR  OBTAINING  IRREDUNDANT 
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LINK( 16384) 

M I NV! 163  841 

CONTINUE 
STARTFU5) 

ENDF       (15) 


r  * 

C    MAIN  PROCEDURE    * 
C  * 

c 

C    READ  PARAMETER  N,  M 


MOS  CELL  CONFIGURATION. 

STORES  THE  LINK  TO  THE  NEXT  VERTEX  WITH  THE  SAME 

WEIRHT  IN  LARGE  CUBE. 

STORES  THE  MINIMUM  VECTOR  OBTAINED  IN  THE  PROCESS 

OF  IMC. 

POINTER  TO  THE  FIRST  VERTEX  WITH  THE  SAME  WEIGHT 
IN  LARGE  CUBE. 

POINTER  TO  THE  LAST  VERTEX  WITH  THE  SAME  WEIGHT 
IN  LARGE  CUBE. 


10  READ(  5,lltEND  =  130  )  N,M 
U  FORMAT!  I2«2XtIl  > 


HEADING 
l»t 


PRINT  12 

12  FORMAT! 
PRINT  13 

13  FORMAT! 
PRINT  14 

14  FORMAT! 
PRINT  15 

15  FORMAT! 
PRINT  16 

16  FORMAT! 
PRINT  17 

17  FORMAT! 
PRINT  18 

18  FORMAT! 
PRINT  19 

19  FORMAT! 


a* t ******************************  ********* 

' t '*  ** 

%•*    DESIGN  OF  IPREDUNDANT  MOS  NETWORK    *• 

i 


i  f  i  a**************************************** 

OS'    X  =  EXTERNAL  VARIABLE'  ) 
•f«    U  =  OUTPUT  OF  MOS  CELL'  I 

Ql  f*  ******************************************       ) 

CHECK  PARAMETER  VALUE 


IF!  M  .FQ.  I  )  GO  TO  21 
IF!  M  .EO.  2  )  GO  TO  22 
IF!  M  .EQ.  3  )  GO  TO  23 
IF!  M  .EO.  4  )  GO  TO  24 
GO  TO  2  5 

21  IF(  N.GE.l  .AND.  N.LE.10  )  GO  TQ  30 
GO  TO  25 

22  IF!  N.GE.l  .AND.  N.LE.9  )  GO  TO  30 
GO  TO  25 

23  IF!  N.GE.l  .ANO.  N.LE.8  )  GO  TO  30 
GO  TO  2  5 

24  IF!  N.GE.l  .AND.  N.LE.7  )  GO  TO  30 

ERROR  IN  PARAMETER  VALUE  i 

25  PRINT    26tN.M 

26  FORMAT!     'QS'INPUT    ERROR    IN    PARAMETER    CARD'tSX.'N    **tI3t'M 
GO    TO    130 


13    ) 
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C    PRINT  PARAMETER  VALUE 
C 

30  PRINT  31, N 
.  31  FORMAT!  •-•t'NUMBER  OF  EXTERNAL  VARIABLES  -',13  ) 
PRINT  32, M 
32  FORMAT(  •  •♦•NUMBER  OF  OUTPUT  FUNCTIONS    -»|I3  ) 
C 

C    INITIALIZE  LABEL!  2**N  )  AND  DCARE!  2**N  I 
C 

NVTEXL  ■  2**N 
DO  50  I  -    1, NVTEXL 
LABEL!  I  )  »  0   «.• 
DCARE!  I  )  =  0  ■ 

50  CONTINUE 

CALL  INPUT!  CS.C130  ) 
C 

C    SET  TIMFR 
C 

CALL  STEPZ!  TIME  ) 

CALL  NCUBE 

11  =  1 

RF  =  16 

TFET  =0 

PRINT  51 

51  FORMAT (  •-•♦•NETWORK  CONFIGURATION') 
60  CALL  CMNL 

IF(  RF  .EQ.  I  )  GO  TO  90 
CALL  CMXL 
C 

C    IF  MMI  (  I  )  =  MXL(  I  I  FOR  EVERY  I  IN  N-CUBE,  CMNLt  CMXL,MPF 
C    CAN  BE  SKIPPED. 
C 

DO  70  I  =  lv NVTEXL 

IF(  MNL(  I  )  .NE.  MXL(  I  I  I  GO  TO  80 
70  CONTINUE 

GO  TO  120 
80  CALL  MPF 

CALL  IMC(  TFET, CIO  I 
11=11+1 
IF(  CS  )  GO  TO  100 
IF(  II  .NE.  RF  )  GO  TO  60 
90  CALL  ASFUN1 

CALL  IMC(  TFET, £10  ) 
GO  TO  125 
100  IF<  RF  -  II  ♦  1  .GT.  M  )  GO  TO  60 
C 

DECIDE  THE  MOS  CELL  CONFIGURATIONS  IN  COMPLETELY  SPECIFIED 
C    FUNCTION  PART 
C 

110  CALL  ASFUN1 

CALL  IMC(  TFET, CIO  ) 
II  =  II  ♦  I 

IF(  II  .LE.  RF  )  GO  TO  110 
GO  TO  125 
C 

C         MNL(     I     )    =    MXL(     I     )    OCCURED    FOR    EVERY    I     IN    N-CUBE 
C 

120  CALL  ASFUN2 

CALL  IMC(  TFET, CIO  ) 

II  =  I  I  ♦  1 

I F I  II  .LE.  RF  J  GO  TO  120 
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MOS 

CELLS 

■  • 

»I3 

) 

MQS 

FETS 

a* 

»I3 

) 

FACTORING 

)' 

) 

!  / 

100. 

ME  »• ,F7.5 

It' 

SEC 

• 

c 

C    PRINT  STATISTICS 
C 

•125  PRINT  126, RF 

126  FORMAT!  •-', 'NUMBER  OF 
PRINT  127,TFET 

127  FORMAT!  •  ', 'NUMBER  OF 
PRINT  129 

129  FORMAT!  •  ', • I  WITHOUT 
CALL  STEPZ(  ITIME  ) 
UTIME  =  (  TIME  -  ITIME 
PRINT  128, UTIME 

128  FORMAT!  '6', 'ELAPSED  TIME  *',F7.2,'  SEC  I 
GO  TO  10 

130  STOP 
END 

£**«******************** 

C  * 

C    SUBROUTINE  NCUBE    * 
C  * 

Q ***********  ************ 

SUBROUTINE  NCUBE 
C    THIS  SUBROUTINE  CONSTRUCTS  NCUBE  ACCOROING  TO  THE  NUMBER  OF  EXTERNA! 
C    VARIABLES.  CHECK  THE  INPUT  VECTOR  IN  BINARY  FORM  ASSIGNED  TO  EACH 
C    VERTEX  ONE  3Y  ONE  AND  LINKS  THE  VERTICES  WITH  THE  SAME  WEIGHT. 

IMPLICIT  INTEGER!  A-Z  ) 

INTEGFR*2   LABEL   ,DCARE   ,MNL     fMXL     t CHAIN   ,FUNC    tLINK 
1  ,MINV 

COMMON       II       ,N        ,M        ,R.F 

1  ,LA3EL<1024)  ,DCARE!1024)  ,MNL!1024>  ,MXL<1024) 

2  tCHAIN(l024)  ,STARTL!U)  tENDLlll)  tFUNC!16384) 

3  ,LINK(16384)  ,MINVU6384)  ,STARTF!151  ,ENDF!15) 
MXVTXL    =    2**N    -     1 

L    =    N    ♦    I 

DO    10    I    =    1,L 
10    STARTL!     I     )    =    0 

I    =    0 
20    X    =    I 

W    =    0 
30    IF!    X.EO.O    I    GO    TO    50 

IF!     X/2*2.EQ.X    )    GO    TO    40 

W    =    W    ♦     1 
40    X    -    X/2 

GO    TO    30 
50    IF!     STARTL!    W    ♦    I    ).EQ.0    I    GO   TO   60 

CHAIN!     ENDL!     W    ♦     I     )     1=1+1 

ENDL!     W    +    1     )     =    I    ♦    1 

GO  TO  70 
60  STARTL!  W  ♦  1  1  =  I  ♦  1 

ENDL!  W  ♦  1  )  >  I  ♦  1 
70  I  =  I  ♦  1 

IF!  I. LE. MXVTXL  )  GO  TO  20 

RETURN 

END 
C************* ********** 

C  *  -  / 

C    SUBROUTINE  INPUT    * 

C  * 

c***************** ****** 

SUBROUTINE  INPUT!  CSt*  ) 
C    THIS  SUBROUTINE  REAOS  DATA  INTO  LABELt  2**N  )  AND  DCARE!  2**N  ) 
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C    ON  INPUT  ERROR,  THE  FOLLOWING  MESSAGE  IS  PRINTEO  OUT  AND  EXECUTION 

C    IS  TERMINATFO. 

C  I  =  FUNCTION  NO.   J  «  CARD  NO.   K  ■  COLUMN  NO. 

IMPLICIT  INTEGER(  A-Z  ) 

INTEGER*2   LABEL   ,DCARE   ,MNL"    »MXL     , CHAIN   ,FUNC    ,LINK 
1  ,MINV 

COMMON       II       ,N        ,M        ,RF 

1  ,LABEL!1024)     ,CCARE!1024)     ,MNL!1024)       ,MXL<1024) 

2  ,CHAINU024>     ,STARTL!ll)       ♦ENDL(ll)         ,FUNC!16384) 

3  ,LINK!16384)     .MINVI 163841      ,STARTF(15)      ,ENDF(15) 
LOGICAL  CS 

INTFGER*2  CHAR(80)   - 

INTFGER*2  BLANK/'  • / ,ZERO/ ' 0 ' / , ONE/ • 1 ' /,0CR/ •* • / 
C    CHARI  80  )        CHARACTER  BUFFER  FOR  ONE  CARD. 
C    CS  IS  SET  IF  THE  GIVEN  FUNCTIONS  ARE  COMPLETELY 

C  SPECIFIED. 

CS  =  .TRUE. 

IF(  2**N/80*80  .EQ.  2**N  )  GO  TO  10 
NCAPD  =  2**N/80  ♦  I 
GO  TO  20 
10  NCAPD  =  2**N/80 
20  00  140  I  =  1,M 

PPINT  30, I 
30     FORMAT!  • 0 « , • FUNCT ION' , I  3  1 
VTEX  =  i 

DO  100  J  =  1,NCARD 
READ  40,CHAR 
40         FORMAT!  80A1  ) 

DC   flO   K   =    1    80 

IF(  CHAR!  K  I  .EQ.  BLANK  )  GO  TO  110 
IF(  CHAR(  K  )  .EQ.  ZERO  I  GO  TO  70 
IF(  CHARI  K  J  .EQ.  ONE  )  GO  TO  50 
IF{  CHAR!  K  I  .EQ.  DCR  )  GO  TO  60 
GO  TO  150 
50  LABEL(  VTEX  )  ■  LABEL!  VTEX  )  ♦  2**(  M  -  I  I 

GO  TO  70 
60  OCAREC  VTEX  )  =  DCARE1  VTEX  I  +  2**(  M  -  I  ) 

CS  =  .FALSE. 
70  VTEX  =  VTEX  ♦  I 

80        CONTINUE 

PRINT  90, CHAR 
90        FORMAT!  »  ,,3X,80A1  I 
100     CONTINUE 

GO  TO  120 
110     IF!  J  .NE.  NCARD  )  GO  TO  150 

IFI  VTEX  .NE.  2**N  ♦  1  )  60  TO  150 
120     PRINT  130, CHAR 
130     FORMAT!  '  ',80Al  ) 
1^»0  CONTINUE 

GO  TO  170 
150  PRINT  160. I , J.K 

160  FORMAT!  •0,tMNPUT  ERROR  IN  DATA  CARD',3X,M 
I  'K  =»,I3  ) 

PFTURN1 
170  OFTURN 
FND 
r *********** *********** 

c  * 

C    SUBROUTINE  CMNL    * 

C  * 

q ********************** 
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SUBROUTINE  CMNL 
C    THIS  SUBROUTINE  IMPLEMENTS  CONDITIONAL  MINIMUM  LABELING, 
C    THIS  SUBROUTINE  CALLS  SUBROUTINE-DSCAN  AND  SUBROUT INE-INCMNT. 

IMPLICIT  INTEGER!  A-Z  I 

INTEGER*2   LABEL   tDCARE   »  MNL     ,MXL     tCHAIN   ,FUNC    ,LINK 
1  ,MINV 

COMMON       II       ,N        ,M        ,RF 

1  ,LAB?L!1024)     ,0CAREI1024)      .MNLI1024)        .MXLU024) 

2  .CHAIN11024)     ,STARTLlll)       tENDL(il)         tFUNC<  16384) 

3  ,LINK(16384)     ,MINV!16384)     ,STARTFC15J      ♦  ENDFUS) 
INTEGER  BWEITJ4) 

C    PTR8  POINTS  TO  THE  VERTEX-B  FOR  WHICH  WE  ARE  SEEKING  FOR  THE 
C    MINIMUM  POSSIBLE  LABEL. 

C    PTRA  POINTS  TO  THE  VERTEX-A  WHICH  HAS  BIGGER  WEIGHT  BY  ONE  THAN 
C    VERTEX-B  AND  CONNECTEO  TO  VERTEX-B  BY  THE  EDGE  IN  N-CUBE. 

FF!  MINIMUM  NUMBER  OF  MOS  CELLS  )  IS  OBTAINED  IN  THE  FIRST  EXECUTION 
C    OF  THIS  SUBROUTINE. 
F  -  M 

MNL  (  STARTL(  N  ♦  1  )  )  =  LABEL!  STARTLt  N  ♦  1  )  ) 
USFFY  =  RF  -  I!  ♦  I 
W  =  N 
10  PTPB  =  STARTL(  W  ) 
20  LBX  =  LABEL!  PTRB  ) 
LBY  =  DCAREl  PTRB  ) 
INC  «  0 
SHIFT  =  1 

CALL  DSCAN{  LBY,USPFY,NDCAREtBWEIT  ) 
DO  50  I  =  l.N 

PTDBX  =  (  PTRB  -  I  )  /  SHIFT 
IF(  PTRBX/2*2  .NE.  PTRBX  )  GO  TO  40 
PTRA  ■  PTRB  ♦  SHIFT 
LAX  =  MNL!  PTRA  ) 
3  0     IF(  LBX  ,GE.  LAX  )  GO  TO  40 
LBX  =  LABEL!  PTRB  ) 

CALL  INCMNT!  I NC» NDCARF , L BX , F»BWE IT  » 
GO  TO  30 
40     SHIFT  =  2  *  SHIFT 
50  CONTINUE 

MNL!  PT"B  )  =  LBX 

TF(  PTRB  ,E0.  ENDL!  W  )  )  GO  TO  60 
PTRB  =  CHAIN!  PTRB  I 
GO  TO  20 
60  W  =  W  -  1 

IF!  W  .GE.  I  )  GO  TO  10 
IF(  I  I  .NE.  1  )  GO  TO  100 
RF  =  0 

MNLX  =  MNL!  1  ) 
70  IF!  MNLX  .EO.  0  I  GO  TO  100 
MNLX  =  MNLX  /  2 
RF  =  RF  ♦  1 
GO  TO  70 
100  RETURN 
END 

Q *********** *********** 

C  * 

C    SUBROUTINE  CMXL    *  "  ' 

C  * 

£********************** 

SUBROUTINE  CMXL 
C    THIS  SUBROUTINE  IMPLEMENTS  CONDITIONAL  MAXIMUM  LABELING. 
C    THIS  SUBROUTINE  CALLS  SUBROUTINE-DSCAN,  SUBROUT I NE-AS IGNL  AND 
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C  SUBROUTINE-DCRMNT. 

IMPLICIT    INTEGER!    A-Z    ) 

INTEGER*2      LABEL       ,DCARE       ,MNL  tMXL  , CHAIN       ,FUNC  ,LINK 

1  ,MINV 

COMMON  II  ,N  ,M  ,RF 

1  ,LABEI<1024)  ,DCARE(1024)  ,MNL(1024)  ,MXL(1024) 

2  ,CHAIN(1024>  tSTARTLflll  ,ENDL<11)  ,FUNC(16384) 

3  ,LINK<16384)  ,MINV<163841  ,STARTF(15)  ,ENDF(15) 
INTEGER     BWEIT14) 

C  PTRB    POINTS    TO    THE    VERTEX-B    FOR    WHICH    WE    ARE    SEEKING    FOR    THE 

C  MAXIMUM    POSSIBLE    LABFL. 

C  PTRA    POINTS    TO    THE    VFRTEX-A    WHICH    HAS    SMALLER    WEIGHT     BY    ONE    THAN 

C         VERTEX-B    AND    CONNECTEDTO    VERTEX-B    BY    THE    EDGE    IN    N-CU8E. 
F    =    M 

LBX    =    LABEL!     STARTLl     I     )     ) 
LBY    =    DCAPEl     STARTLl     1)1 
USPFY    =    RF    -    II     +1 
CALL    ASIGNIC    U SPFY, F ,LRX,LBY    I 
MXLC     STARTLt     I     )     )     =    LBX 
L    =    N    ♦    1 
DO    60    W    =    2,L 

PTRR    =    STARTLl     W     ) 
20  LbY    =    DCAREI     PTRB     ) 

CALL    DSCAN(     LB Y, USPFY ,NDC ARE ,BWE I T    ) 

LBX    =    LABEL(     PTRB     ) 

LBY    =    DCAREI     PTRB     ) 

CALL    ASIGNK    USPFY,  F  ,  LB  X, LBY    ) 

LBXX  =  LBX 

INC  =  0 

SHIFT  =     1 

DO  50  I  =  1,N 

PTRBX  =  (  PTRB  -11/  SHIFT 
IF(  PTPBX/2*2  .EQ.  PTRBX  )  GO  TO  40 
PTRA  =  PTRB  -  SHIFT 
LAX  =  MXL(  PTRA  ) 
30        IFI  LBXX  .LE.  LAX  )  GO  TO  40 
LRXX=LBX 

CALL  DCRMNT<  I NC, NDC ARE, L BXX, F,BWE  IT  ) 
GO  TO  30 
40        SHIFT  =  2  *  SHIFT 
50     CONTINUE 

MXL(  PTRB  )  =  LRXX 

IF<  PTRB  .FQ.  ENDL(  W  )  )  GO  TO  60 
PTRP  =  CHAIN(  PTRB  ) 
r.n  TO  20 
60  CONTINUE 
RF  TURN 
END 

£*********** ************ 

C  * 

C    SUBROUTINE  DSCAN    * 

c  * 

c*********** ************ 

SURPOUTINF  DSCAN(  LBY.USPF Y, NDC ARE , BWE I T  ) 
C    THIS  SUBROUTINE  SCAN  THE  OCARE  WHICH  IS  ASSIGNED  TO  EACH  VERTEX 

AND  OBTATN  THE  NUMBER  OF  DONTCARE  BITS  AND  THE  WEIGHT  ASSIGNED 
C    TO  rACH  DONTCARE  BIT. 

MPl  IC  IT  INTEGER!  A-Z  ) 
riTtf.FR  BwEITC.  ) 
C    pwcIT(  4  )        STHRES  THE  WEIGHT  ASSIGNED  TO  EACH  DONTCARE  BIT. 
Nf'CARF  STORFS  THE  NUMBER  OF  DONTCARE  BITS  IN  ONE  VERTEX. 
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NDCARE  -  0 

1  =  1 

SHIFT  ~    1 

DO  20  I  *  1»USPFY 

IF(    LBY    „EQ.    0    )    GO   TO    30 
IF(    L8Y/2*2    .EQ.    LBY    >    GO   TO    10 
NDCARE    =    NDCARE    ♦    I 
BWEIT(    NDCARE    )    *    SHIFT 
10  SHIFT    =    2    *    SHIFT 

LRY  =  LBY  /  2 
20  CONTINUE 
30  RETURN 
END 

£******?**«************** 

C  * 

C    SUBROUTINE  INCMNT    * 
C  * 

£*****#****************** 

SUBROUTINE  INCMNTI  I NC .NDCARE, LBXtMt BWEIT  I 
C    THIS  SUBROUTINE  INCREMENTS  THE  LABEL  ASSIGNED  TO  THE  VERTEX  POINTED 
C    BY  POINTER-PTRB. 

IMPLICIT  INTEGER!  A-Z  I 

INTEGER  BWEITl  4  \ 
C    INC  IS  THE  COUNTER  WHICH  HAS  SPECIFIED  WEIGHT  FOR  EACH  BIT. 

INC  =  INC  M 

INCX  =  INC 

J  =  1 

SHIFT  =  2  **  (  M  -  NDCARE  ) 
10  IF (  INCX  .EQ.  0  )  GO  TO  40 

IFl  INCX/2*2  .EQ.  INCX  )  GO  TO  30 

IF(  J  .IE.  NDCARE  )  GO  TO  20 

LBX  =  LBX  ♦  SHIFT 

GO  TO  30 
20  LBX  =  LBX  ♦•  BWEIT(  J  I 
30  J  =  J  ♦  1 

SHIFT  =  2  *  SHIFT 

INCX  =  INCX  /  2 

GO  TO  10 
40  RETURN 

END 
C *************  *********** 

C  * 

C    SUBROUTINE  ASIGN1    * 
C  * 

C************* *********** 

SUBROUTINE  ASIGNK  USPFY.M.LBX, LBY  ) 
C    THIS  SUBROUTINE  ASSIGNS  ONE  TO  THE  EVERY  UNSPECIFIED  BIT  OF  LABELS 
C    IN  N-CUBE. 

IMPLICIT  INTEGER!  A-Z  ) 
SHIFT  =  I 

IF(  USPFY.GT.M  )  GO  TO  30 
DO  20  I  =  l.USPFY 

IF(  LBY/2*2  .EQ.  LBY  )  GO  TO  10 
LBX  =  LBX  ♦  SHIFT 
10     SHIFT  =  2  *  SHIFT 

LBY  =  LBY  /  2  -  ' 

20  CONTINUE 
GO  TO  70 
30  00  50  I  =  1,M 

IFl  LBY/2*2  .EQ.  LBY  I  GO  TO  40 
LBX  =  LBX  ♦  SHIFT 
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40     SHIFT  »  2  *  SHIFT 
LBY  =  LBY  /2 

50  CONTINUE 
I  =  M  +  I 
DO  60  I  -    UUSPFY 

LBX  =  L8X  ♦  SHIFT 
SHIFT  =  2  *  SHIFT 
60  CONTINUE 
70  RETURN 
END 

£************************ 

c  *  "  , 

C    SUBROUTINE  DCRMNT    * 

C  * 

f-*****:******************* 

SUBROUTINE  DCRMNT {  I NC t NOCARE, LBXX ,M, BWE IT  ) 
C    THIS  SUBROUTINE  DECREMENT  THE  LABEL  ASSIGNED  TO  THE  VERTEX  POINTED 
C    BY  POINTER-PTRB 

IMPLICIT  INTEGER(  A-Z  ) 

INTEGER  BWEITI  4  I 

INC  =  INC  ♦•  1 

INCX  =  INC 

J  =  1 

SHIFT  =  2  **  (  M  -  NDCARE  ) 
10  IF  f  INCX.EO.O  )  GO  TO  40 

I*<  INCX/2*2  .EO.  INCX  )  GO  TO  30 

IF(  J. LE. NDCARE  )  GO  TO  20 

L^XX  =  LBXX  -  SHIFT 

GO  TO  30 
20  LBXX  =  LBXX  -  BWEITI  J  ) 
30  J  =  J  ♦  I 

SHIFT  =  2  *  SHIFT 

INCX  =  INCX  /  2 

GO  TO  10 
40  PFTURN 

END 
C ******* ************** 

C  * 

C    SUBROUTINE  MPF    * 
C  * 

Q ************* ******** 

SUPPOUTINE  MPF 
C    THIS  SUBROUYINE  OBTAINS  MAXIMUM  PERMISSIBLE  FUNCTION  FROM  MINIMUM 
C    LABEL  AND  MAXIMUM  LABEL  ANO  STORES  THE  RESULT  TO  FUNCI  NVTEXF  ). 

IMPL  If.  IT  INTEGER!  A-Z  ) 

INTEGER*2   LABEL   ,DCARE   ,MNL     fMXL     »CHAIN   ,FUNC    ,LINK 
1  ,MINV 

COMMON      II      (N       fM       ,RF 

1  ,LABEL(1024)     ,DCARE(10241     .MNH1024)        ,MXLI1024) 

2  ,CHA[N<1024)      .STAPTL(ll)       tENDL(ll)         ,FUNC(16384) 

3  tLINK(l6384)      »MINV(16334)      ,STARTF(15)       ,ENDF(15) 
MXVTXL  =  2**N  -  I 

NVTEXF  =  2**1     N  ♦  II  -  1  ) 

PSHIFT  =  2**(  RF  -  II  J 

LSHJFT  =  2**(  II  -  1  ) 

DO  10  I  =  1, NVTEXF  "  / 

10  FUNCt  I  )  =  0 

I  =  0 
20  mm  x  ■  MNLl  I  ♦  I  »  /  RSHIFT 

MXLX  =  MXL(  1*1)/  RSHIFT 

J  =  I*LSHIFT  ♦  MNLX/2 
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.EQ.MNLX    I    GO    TO   30 
.EO.    MXLX    )    GO    TO    40 
»    2 

•NE.    MXLX    )    GO    TO    40 
*    1 


C* 

c 
c 
c 
c* 

c 
c 
c 
c 

c 
c 

r 

c 

c 
c 
c 
c 
c 
r 
c 
c 
c 
c 


IF(     MNLX/2*2 

IF<     MXLX/2*2 

FUNC(     J    ♦    1     ) 

GO    TO    40 
30    IF(     MXLX/2*2 

FUNC(    J    ♦    i     ) 
40    I    ■    I    ♦    1 

IF|     I.LE.MXVTXL     )    GO    TO    20 

RCTURN 

END 
***************** 

SUBRPUTINE  IMC 

***************** 
SUBROUTINE  IM 
THIS  SUBROUTINE 

THE  MAXIMUM  PER 
THIS  SUBROUTINE 


STEP-I 
STEP-2 


STFP-3 
STEP-4 

STEP-5 

STEP-6 


CONS 

ASSI 

THAT 

F<  Y 

08TA 

OBTA 

COVE 

OBTA 

IN  S 

STOR 

THIS 

IF  F 

IN  P 

AIRE 

THE 

IMPLICIT  INTE 

INTEGER*2   LA 

1  ,MI 

COMMON       II 

1  ,LAPEL(1024 

2  ,CHAIN(1024 

3  .LINKI16384 
INTEGER  X(  10 

I 
2 
3 

INTEGER 
1 
2 
3 

INTEGER 

LOGICAL 

LOGICAL 
PRTAPY(  20 


*** 

* 

* 

* 

*** 

C(  TFET.*  ) 

OBTAINS  IRRE 
MISSIBLE  FUNC 

IS  MADE  OF  T 
TRUCTS  LARGC- 
GNS  "0"  OR  "1 

THERE  EXISTS 

)  =  "1". 
INS  THE  SFT  0 
INS  THE  SUBSE 
RS  EVERY  VERT 
INS  AN  IRREOU 
TEP-4. 
ES  THE  RESULT 

STEP  CONTAIN 
UNCTION  VALUE 
REVICUS  STEPS 
ADY  STORED  IN 
CONTENTS  OF  A 
GER(  A-Z  J 


DUNOANT  MOS  CELL  CONFIGURATION  FROM 
TION. 

HE  FOLLOWING  SIX  STEPS. 
CUBE  FUNC(  2**(N  ♦  II  -  1)  ). 
"  TO  THE  VERTEX-X  WITH  OONTCARE  SUCH 
NO  VERTEX-Y  SATISFYING  Y  >  X  AND 

F  MINIMUM  VECTORS. 

T    OF    THE    SET    OF    MINIMUM    VECTORS    WHICH 

EX  WITH  ORIGINAL  "0". 

NDANT  SUBSET  FROM  THE  SUBSET  OBTAINED 

TO  THE  VERTEX  IN  N-CUBE  (  LABEL  ). 

S  ERROR  CHECKING  ROUTINE. 

S  FRCM  MOS  CELL  CONFIGURATION  OBTAINED 
ARE  DIFFERENT  FROM  THE  FUNCTION  VALUES 
ARRAY-LABEL.  ERROR  MESSAGE  AND 

RRAY-LABEL  IS  PRINTED. 


BEL 
NV 


DCARF   rMNL 


MXL 


♦CHAIN   ,FUNC 


LINK 


)  / 


U<  10  )  / 


iN 

TDCA<?E{  I 
.STARTLC 
f  MINV( 16 
i  _  • 


.M 

024) 

11) 

384) 

—  X2'  , 

— X7» ,• — X8' 

— U2*  ,  • — U3« 


,RF 

,MNL(1024) 
,ENDL{  11) 
,STARTF( 15) 
•  —  X3'  .  • —  X4» 


,MXL(10241 

,FUNC(16384) 
.ENDFU5) 
. — x5' 


•  — i 


BLANK/ 

F0RM3/  • — |  • 
F0RM6/  '-0   • 

F0P.M9/  «0 • 

PRTARYl  20  ) 
ERROR 
REPEA 
) 


, U4»,' — U5 

i — U7# t ' — US' .  • — U9»,«-U10 
*  •  /  ,  F0RM2/ 

•  /  ,  F0RM5/ 

•  /  ,  F0RM8/ 


F0RM1/ 
F0RM4/  •    | 
F0RM7/  • 


=% 


PRINT  BUFFER  FOR  PRINTING  OBTAINED  MOS  GATE 

CONFIGURATION. 

ERROR  IS  SET  IF  ERROR  OCCURS  IN  CONSTRUCTING  NETWORK. 

INPD  =  N  «■  I  I  -  1 
NVTEXL  =  2  **  N 
MXVTXL  =  NVTEXL  -  I 
NVTEXT  =  2  **  INPD 
MXVTXF  =  NVTEXF  -  1 
LSHIFT  =  ?  **  (  II  -  1  ) 
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RSHIFT  «  2  **  (  RF  -  II  ) 

ERROR  =  .FALSE. 


c 

c 

1    STE 

c 

L    ■    INPD    ♦    1 
DO    10    I    =    1,L 

10 

STARTFI     I    )    =    0 
I    =    0 

20 

Y    *    I 

W    =    0 

30 

IF<     Y    .EQ.    0    )    GO    TO  -50 

Ic{    Y/2*2    -EQ.    Y    )    GO    TO   40 

W    =    W    «-     1 

40 

Y    =    Y    /    2 

GO    TO    30 

50 

IF(     STARTF<     W    ♦    I    ).E0.0    1    GO    TO   60 
LINKC     FNDF(     W    +    1     )     »     *    I    ♦     I 
ENDFl     W    ♦    1     1    »     I     ♦    1 

GO    TO    70 

60 

STARTF(     W    «■    1     )    «    I     ♦    1 
ENDFl     W    ♦    1     )    =     I     ♦    1 

70  I  =  I  «■  1 

IF(  I.LE.MXVTXF  )  GO  TO  20 

STEP-2 

IF(  FUNCI  STARTF(  INPD  ♦  I  )  )  ,NE.  0  1  GO  TO  110 
FUNCt  STARTF(  INPD  ♦  1  )  )  =  3 
110  W  =  INPD 
120  PTRB  =  STARTF(  U  ) 

130  IFI  FUNC(  PTRB  J.NE.O  I  GO  TO  160 
SHIFT  =  1 
DO  150  I  =  1,INPD 

PTRBX  =  (  PTR8  -  1  )  /  SHIFT 
IF[  PTRBX/2*2  .NE.  PTRBX  )  GO  TO  140 
PTRA  =  PTRB  +  SHIFT 
IF(  FUNC(  PTRA  J.NE.2  )  GO  TO  140 
FUNC(  PTRB  )  =  2 
GO  TO  160 
140     SHIFT  =  2  *  SHIFT 
150  CONTINUE 

FUNC(  PTRB  I  =  3 
160  IF(  PTRB.EO.ENDF(  W  )  )  GO  TO  170 
PTRB  -    LINK(  PTRB  ) 
GO  TO  130 
170  W  »  W  -  1 

IF(  W.GE.l  )  GO  TO  120 

STEP-3 

K  =  I 

NMINV  =  0 
W  =  INPD  ♦  1 
210  PTRR  =  STARTFI  W  ) 
220  IF(  FUNCf  PTRB  J.EQ.2 
SHIFT  =  I 
00  ?40  I  =  1,INPD 

PTPBX  =  I  PTRB  -  I 
ic(  PTRBX/2*2  .EQ. 
PTRA  =  PTRB  -  SHIFT 


)  GO  TO'250 


)  /  SHIFT 

PTRBX  I  GO  TO  230 
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IF(  FUNCI  PTRA  ).NE.2  )  GO  TO  250 
230     SHIFT  »  2  *  SHIFT 
240  CONTINUE 

'    MINV(  K  J  »  PTRB  -  1 
NMINV  »  NMINV  ♦  1 
K  =  K  ♦  1 
250  IF(  PTRB.EQ.ENCF(  H  )  )  GO  TO  260 
PTRB  =  LINKl  PTRB  J 
GO  TO  220 
260  W  =  W  -  1 

IF(  W.GE.l  I  GO  TO  210 
C 

C    STEP-4 
C 

DO  310  I  =  1,NVTEXF 
310  I INK(  I  )  *  0 

I  =  1 
320  REPEAT  =  .FALSE. 

DO  ?40  J  =  i,NVTEXF 

IF(  FUNCI  J  I.NE.l  )  GO  TO  340 
JX  =  J  -  1 
MINX  =  MINVI  I  ) 
CALL  COMPRl  JX,MINX,C330  I 
LINK!  J  )  =  LINK!  J  )  ♦  1 
GO  TO  340 
330     IF(  LINK!  J  J.NE.O  )  GO  TO  340 

REPEAT  =  .TRUE. 
340  CONTINUE 

IF!  REPEAT  )  GO  TO  350 
NSIRR  =  I 
GO  TO  360 
350  I  »  I  ♦  1 
GO  TO  320 
360  CONTINUE 
C 

C    STEP-5 
C 

L  =  NSIRR  -  1 
DO  450  I  =  1,L 

DO  410  J  =  l,NVTEXF 

IF(  FUNCI  J  J.NE.l  I  GO  TO  410 
JX  =  J  -  1 
MINX  ■  MINVI  I  J 
CALL  COMPRl  JX, MINX, 6410  ) 
LINK!  J  )  =  LINK!  J  )  -  1 
IF(  LINK!  J  J.EO.O  I  GO  TO  420 
410     CONTINUE 

GO  TO  440 
420     IF(  FUNCI  J  l.NE.l  I  GO  TO  430 
JX  =  J  -  1 
MINX  =  MINVI  I  ) 
CALL  COMPRl  JX, MINX, £430  ) 
LINK(  J  )  =  LINK!  J  »  ♦  1 
430     J  =  J  -  1 

IF(  J.CE.l  )  GO  TO  420 
GO  TC  450 
44  0     MINVI  I  )  =  -1 
450  CONTINUE 
C 

C    PRINT  MOS  GATE  CONFIGURATION 
C 
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NIRR 
MXFE 
00  4 


460 

470 
480 

490 
tlOO 


NFE 


CONT 

p?  1  n 

FORM 
I  = 
L  = 
DO  4 
I 


=  0 

T  = 
90  I 
F(  M 
i  RR 
INX 

-  0 
F(  M 
F(  M 
FET 
INX 
0  TO 
FET 
Fi  N 
Xi:ET 
INUE 
T  41 
ATI 
0 

?  * 
290 
F(  L 


=  ltNSIRR 
1NVI  I  J  .LT. 
=  NIRR  *    I 
*    MINVI  I  ) 


0  )  GO  TO  490 


INX  .EQ. 
INX/2*2  , 

«  NFET  ♦ 
=  MINX  / 

460 
■  TFET  ♦ 
FET  .LE. 

=  NFET 


0  )  GO  TO  480 

EO.  MINX  J  GO  TO  470 

1 

2 


NFET 
MXFET 


I  GO  TO  490 


00 
•0' 


J 


NIRR  -  ?. 
LINE  =  1,1 
INE/2*2  .EO. 


LINE  J  GO  TO  4240 


PRINT  000  LINE 


4101 


4110 


4111 
4120 


4130 
4140 


4150 
4160 
4170 


4180 
4190 


I  = 
IF( 
I  F  { 

PRT 
GO 
IF  I  NI 
PRTARY 
GO  TO 
PRTARY 
K  = 
SHI 
DO 


CON 
IF( 

Nl 
N2 
00 


CON 

K  = 

IF( 

PPT 

GO 

IFI 

PMT 

PRT 
MX3 


I  ♦ 
MIN 
LIN 
ARY( 
TO  4 
RR  . 
(  1 
4120 
(  1 

0 
FT  = 
4140 
MINX 
IF( 
K  = 
PRTA 
SHIF 
TINU 

I  I 
=  N 
=  N 
A  160 
MINX 
IF  ( 
K  = 
PRTA 
SHIF 
TINU 
K  ♦ 
K  . 
ARY( 
TO  4 
L  IN 
«JY( 
AHY( 
■  M 


I 

VI  I 
E  .EQ 

I  I 
120 
EQ.  1 
)  =  F 


)  .LT.  0  )  GO  TO  4101 
.  NIRR  )  GO  TO  4110 
=  F0RM1 

I  GO  TO  4111 

0RM2 


)  =  F0RM9 


2  ** 
J  = 

=  MI 
MINX/ 
K  ♦  1 
RYl  K 
T  =  S 
E 
.EO. 

♦  1 

♦  I  I 
J  = 
=  MI 

MINX/ 
K  ♦  1 
RYl  K 
T  =  S 
E 

1 
GT.  M 

K  ♦ 
170 
F  .EQ 

MXFt 

MXFE 
XFFT 


I  N  +  II  -  2  ) 

ItN 

NVI  I  )  /  SHIFT 

2*2  .EQ.  MINX  )  GO  TO  4130 


♦  I  )  =  X( 

HIFT  /  2 


J  ) 


1  I  GO  TO  4170 

-  1 

Nl  ,N2 

NVI  I  )  /  SHIFT 

2*2  .EQ.  MINX  1  GO  TO  4150 

♦  I  )  ■  UN  -  N  1 
HIFT  /  2 


XFET  )  GO  TO  4180 
1  )  =  F0RM7 

.  NIRR  )  GO  TO  4210 
T  ♦  2  )  =  F0RM3 
T  +  3  1  =  BLANK 
♦  3 
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PRINT  4200,!  PRTARY!  JJ  |,JJ  ■  i,MX3  ) 
4200     FORMAT!  •  »,12X,20A4  I 
GO  TC  4290 

4210  IF(  NIRR  .EO.  I  I  GO  TO  4211 
PRTARY(  MXFET  ♦  2  )  =  F0RM5 
PRTARY!  MX^ET  ♦  3  )  *  F0RM6 
GO  TO  4220 

421 1  PRTARY!  MXFET  ♦  2  »  »  F0RM7 
PRTARY!  MXFET  ♦  3  I  «  F0RM6 

4220     MX3  =  MXFET  ♦  3 

PRINT  4230,11, (  PRTARY!  JJ  I ,JJ  »  l,MX3  I 
4230  FORMAT!  •  ','MOS  CELL' , 1 2, 2X ,20A4  I 

GO  TO  4290  ! 

C 

C    PRINT  EVEN  LINE 
C 
4240     IF!  LINE  .EO.  NIRR  )  GO  TO  4250 
PRTARY!  1  1  =  F0RM1 
GO  TC  4260 
4250     PRTARY!  1  )  =  F0RM2 
4260     DO  4270  J  =  1, MXFET 
427C     PRTARY!  J  +  1  )  =  BLANK 

IF!  LINE  .EQ.  NIRR  »  GO  TO  4280 
PRTARY!  MXFET  +21=  F0RM4 
PRTARY!  MXFET  ♦  3  J  *  BLANK 
GO  TO  4190 
4230     PRTAPY!  MXFET  ♦  2  )  =  F0RM8 
PPTARY!  MXFET  +  31=  F0RM6 
GO  TO  4220 
4290  CONTINUE 

IF!  NIRR  .EQ.  NMINV  )  GO  TO  610 
C 

C    STEP-6 
C 

C    THE  SIZE  OF  IRRFDUNDANT  SUBSET  IS  NOT  EQUAL  TO  THE  SIZE  OF 
C    MINIMUM  VECTOR  SET. 


C 


I  =  0 
510  J  =  I*LSHIFT  ♦  LABEL!  I  ♦  1  )/(  2  *  RSHIFT  I 

DO  520  K  ■  1,NSIPR 

IF!  MINVl  K  )  .LT.  0  I  GO  TO  520 

JX  *  J 

MINX  =  MINVl  K  » 

CALL  COMPR!  JX,MI NX, 652 0  ) 

GO  TO  540 
520  CONTINUE 

IF!  RF  -  II  ♦  1  .GT.M  »  GO  TO  530 
DCARX  =  DCARE!  1*11/  RSHIFT 
IF!  DCARX/2*2  .NE.  DCARX  }  GO  TO  530 
LABELX  =  LABEL!  I  *■  1  I  /  RSHIFT 
IF{  LABELX/2*2  .NE.  LABELX  )  GO  TO  550 
ERROR  =  .TPUE. 
GO  TO  550 
530  LABEL!  I  ♦  1  )  »  LABEL!  I  ♦  1  )  ♦  RSHIFT 

GO  TO  550 
540  IF!  PF  -  II  ♦  I  .GT.  M  )  GO  TO  550 
OCAPX  =  DCARE!  I  ♦  1  )  /  RSHIFT 
IF!  DCARX/2*2  .NE.  DCARX  )  GO  TO  550 
LABFLX  =  LABEL!  1*11/  RSHIFT 
IF!  LABELX/2*2  .EQ.  LABELX  )  GO  TO  550 
ERROR  =  .TRUE. 
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550  1=1*1 

IF(  I  .LE.  MXVTXL  )  GO  TO  510 

GO  TO  660 
C 

C    THE  SIZE  OF  IRREDUNDANT  SUBSET  IS  EQUAL  TO  THE  SUE  OF 
C    MINIMUM  VECTOR  SET. 
C 

61C  I  =  0 

620  J  =  I*LSHIFT  ♦  LABEL!  14-1  l/I  2  *  RSHIFT  I 

IFC  FUNC(  J  ♦  I  )  .EO.  2  )  GO  TO  630 

IF{  RF  -  II  +  1  .GT.  M  )  GO  TO  650 

DCARX    =    DCAREI     1+11/    RSHIFT 

IF(     DCARX/2*2    -NE.    DCARX    )     GO    TO   650  ■ 

LABELX    =    LABEL(     I     +     I     )     /    RSHIFT 

I F (     LABELX/2*2    .EQ.    LABELX     »    GO    TO   650 

ERROR    =    .TRUE. 

GO    TO    650 
630     IF(     RF    -     II     +1     .GT.    M    )     GO    TO    640 

OCAPX    =    OCAPFl     I     +    I     )     /    RSHIFT 

IM    0CARX/2*2    .NE.    DCARX    )     GO   TO    640 

L/'BFLX    =    LABEL(     I     ♦    I     )     /    RSHIFT 

IF(     LA8ELX/2*2    .NE.    LABELX    )    GO    TO    650 

ERROR    =    .TRUE. 

GO  TD  650 
640  LABEL!  I  +  I  )  =  LA8EL!  I  +  1  )  +  RSHIFT 
650  I  =  I  +  1 

IF  I  I  .LE.  MXVTXL  i    GO  TO  620 
660  IF(  ERROR  )  GO  TO  670 

GO  to  690 
670  P«  TNT  68C,  II 
680  FOR^ATI  '0', 'ERROR  IN  CONSTRUCTING  MOS  GATE*.I3  ) 

CALL  PRINT!  LABEL, NVTEXL  ) 

RETURNl 
690  RETURN 

END 

£*********************** 

C  * 

C    SUBROUTINE  COMPR    * 
C  * 

Q** ************** ******* 

SUBROUTINE  COMPR(  JX.MINX,*  ) 
C    THIS  SUBROUTINE  COMPAIRS  THE  SIZE  OF  TWO  VECTORS  JX  AND  MINX 
10  IF!  MlNX/2*2  .EO.  MINX  )  GO  TO  20 

IF!  JX/2*2  .EO.  JX  I  RETURNl 
20  JX  =  JX  /  2 

MINX  =  MINX  /  2 

IF!  JX.NE.O  .OR.  MINX.NE.O  I  GO  TO  10 

RETURN 

END 
c *********** ************* 

c  * 

C    SUPROUTINE  ASFUN1    * 
C  * 

c ***************** ******* 

SUBPOUTINE  ASFUN1 
C    THIr,  SUBPQUTINE  ASSIGNS  THE  LABEL  OR  FUNCTION  VALUE  FROM  N-CUBE 
C    TO  LAPGF-CUBE!  FUNC  ). 

I  MPl  If,  I  T  INTEGER!  A-Z  ) 

I\TFGER*2   LABEL   .DCARE   ,MNL     ,MXL     .CHAIN   ,FUNC    .LINK 
I  ,MINV 

COMMON      II       , N       ,M       ,RF 
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1  ,LABELfl024)     ,DCAREU024|     ,MNL(1024)        fMXLC1024) 

2  ,CHAIN<1024)     .STARTLlll)      tENDLlll  ,FUNC<16384) 

3  .LINK116384)     ,MINV(16384)      tSTARTF(15)      ,ENDFU5) 

MXVTXL  =  2**N-1 

NVTEXF  =  2**(  N  ♦  II  -  I  1 

RSHIFT  =  2**1  RF  -  II  I 

LSHIFT  =  2**1  II  -  I  ) 

DO  10  I  =  I, NVTEXF 
10  FUNCl  I  )  =  0 

I  =  0 
20  LABELX  =  LABEL!  I  ♦  I  I  /  RSHIFT 

DCARX  ■  DCAREl  I  ♦  i  )  /  RSHIFT 

J  =  I*LSHIFT  ♦  LABELX/2  ' 

IFl  DCARX/2*2  .NE.  DCARX  )  GO  TO  40 

IFl  LABELX/2*2  .EO.  LABELX  »  GO  TO  30 

FUNCl  J  ♦  1  )  =  2 

GO  TO  40 
30  FUNC(  J  ♦  1  I  »  1 
40  I  *  I  «■  1 

I F (  I. LE. MXVTXL  )  GO  TO  20 

RETURN 

END 

c  * 

C    SUBROUTINE  ASFUN2    * 
C  * 

£*******£**************** 

SUBPOUTINE  ASFUN2 
C    THIS  SUBROUTINE  ASSIGNS  THE  LABEL  OR  FUNCTION  VALUE  FROM 
C    MNL(  N-CUBE  )  TO  LARGE-CUBE(  FUNC  I. 

IMPLICIT  INTEGER!  A-Z  I 

INTEGEP.*2   LABEL   »OCARE   ,MNL     tMXL     tCHAIN   tFUNC    ,LINK 
1  ,MINV 

COMMGN       II       ,N        iM        ,RF 

1  ,LABELU024)     ,DCARE(1024)      ,MNL<1024>        ,PXL(1024) 

2  .CHAINU024)     ♦STARTL(ll)       ,ENDL<11)         rFUNC<16384) 

3  .LINKU6384)     ,MINV(16384)      »STARTF(15)      tENDF<15) 
MXVTXL  =  2**N  -  I 

NVTEXF  =  2**1  N  ♦  I  I  —  1  I 

RSHIFT  =  2**<  RF  -  II  ) 

LSHIFT  =  2**(  II  -  1  ) 

DO  10  I  =  1, NVTEXF 
10  FUNC(  I  )  =  0 

I  =  0 
20  MNLX  =  MNL(  I  ♦  I  )  /  RSHIFT 

J  =  I*LSHIFT  ♦  MNLX/2 

IF(  MNLX/2*2  .EQ.  MNLX  )  GO  TO  30 

FUNCl  J  ♦  1  I  =  2 

GO  TO  40 
30  FUNCl  J  ♦  1  )  =  1 
40  I  =  I  ♦  i 

IF(  I. LE. MXVTXL  )  GO  TO  20 

PETURN 

END 
£***** ****** ************ 

c  *  -  / 

C    SUBROUTINE  PRINT    * 

C  * 

q***** ****** ************ 

C    THIS  SUBROUTINE  PRINTS  THE  CONTENTS  OF  AN  ARRAY!  LABELt DCARE , ETC.  I. 

C    ARAY  IS  THE  ARRAY  TO  BE  PRINTED. 
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THE  NUMBER  OF  THE  ELEMENTS  IN  THE  ARRAY  TO  BE 


C    SIZE  INDICATES 
C    PRINTFO. 

SUBROUTINE  PRINT!  ARAY.SIZE  ) 
■     IMPLICIT  INTEGER*  A-Z  1 

INTEGER*2  ARAY( 1 0241 t PBUFC 16 ) 
PRINT  10 

FORMAT!  «0« , »*****DUMP*****» 
00  60  I  =  I, SIZE 
ARAYX  =  ARAY(  I  I 
J  =  16 

IF(  AKAYX/2*2  .EQ. 
PBUF(  J  )  =  1 
40 

J  )  =  0 
=  ASAYX  /  2 
-  I 

GE.l  )  GO  TO  20 
50,(  PBUF(  JJ  J,JJ  *  1,16  ) 
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20 


30 
40 


ARAYX  )  GO  TO  30 


50 
60 


GO  TO 

PRUF  ( 

ARAYX 

J  =  J 

IF(  J, 

PRINT 

format! 

CONTINUE 

RETURN 

END 
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